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ABSTRACT 

The giant southern stream (CSS) is the most prominent tidal debris feature in M31's steUar halo 
and covers a significant fraction of its southern quadrant. The GSS is a complex structure composed 
of a relatively metal-rich, high surface-brightness "core" and a lower metallicity, lower surface bright- 
ness "envelope." We present spectroscopy of red giant stars in six fields in the vicinity of M31's GSS 
(including four new fields and improved spectroscopic reductions for two previously published fields) 
and one field on Stream C, an arc-like feature seen in star-count maps on M31's southeast minor axis 
at i? ~ 60 kpc. These data are part of our on-going SPLASH survey of M31 using the DEIMOS instru- 
ment on the Keck II 10-m telescope. Several GSS-related findings and measurements are presented 
here. We present the innermost kinematical detection of the GSS core to date (i? = 17 kpc). This 
field also contains the inner continuation of a second kinematically cold component that was originally 
seen in a GSS core field at i? ^ 21 kpc. The velocity gradients of the GSS and the second component 
in the combined data set are parallel over a range of AR — 7 kpc, suggesting that this may represent 
a bifurcation in the line-of-sight velocities of GSS stars. We present the first kinematical detection of 
substructure in the GSS envelope (S quadrant, i? ~ 58 kpc). Using kinematically identified samples, 
we show that the envelope debris has a ~ 0.7 dex lower mean photometric metallicity and possibly 
higher intrinsic velocity dispersion than the GSS core. The GSS is also identified in the field of the 
M31 dwarf spheroidal satellite And I; the GSS in this field has a metallicity distribution identical 
to that of the GSS core. We confirm the previous finding of two kinematically cold components in 
Stream C, and measure intrinsic velocity dispersions of ~ 10 and ~ 4 km s~^. This compilation of the 
kinematical (mean velocity, intrinsic velocity dispersion) and chemical properties of stars in the GSS 
core and envelope, coupled with published surface brightness measurements and wide-area star-count 
maps, should improve constraints on the orbit and internal structure of the dwarf satellite progenitor. 
Subject headings: galaxies: halo — galaxies: individual (M31) — stars: kinematics — techniques: 
spectroscopic 



1. INTRODUCTION 

Galaxy merg ers play a key role in galaxy formation and 
evolution (e.g., 'Scarle fc Zi nnlfTOTl] White fc Reed flOTl 
[Abraham et al.lil996: ,Spring'el et al.ll2005h . Observations 
have uncovered abundant evidence of the ubiquity of mi- 
nor mergers in the form of tidal streams. Those dis- 
covered in ou r own Milky Way (MW ) include the Magel- 
lanic stream (Mathewson ct al."1974'), Mono ceros stream 
(jYannv et al.,.2003: Roc ha-Pinto et aLll2003D, and Sagit- 
tarius Stream JIbata et al. 1994t iMajewski et all [20031 : 
iNewberg et al.l I2003D . Detailed studies of tidal debris 
features can provide insight into statistical models of hi- 
erarchical structure formation, the structure and merger 
history of an individual galaxy, and the properties of 
dwarf satellite systems. Due to their sparse stellar den- 
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sity and long dynamical times, stellar halos are ideal for 
investigating the merger history of an individual galaxy 
in detail as it is possible for higher density tidal debris 
features to remain identifiable in phase-space for billions 
of years. 

Observations of the Galaxy and theoretical arguments 
indicate that present-day dwarf satellite galaxies are dif- 
ferent from the earlier generations of dwarf galaxies that 
merged to form the bulk of the virialized stellar halo. 
Stars in the classical satellites of the Galaxy on aver- 
age have lower relative abundances of alpha elements 
than stars of similar metallicities in t he Galaxy's halo 
(lFuhrmannllT998HShetrone et al.| [200ll. [20031 : iVenn et all 
120041: iGeisler et al.l I2007D . although the relative abun- 
dances of alpha elements are similar in the lowest metal- 
ficity stars (Kirby et al. 2009; Frcbel et al. 2009). This 
has been shown to be a natural consequence of the hierar- 
chical formation framework: simulations of stellar halo 
formation indicate that the majority of stellar mass in 
the halo was contributed by massive satellite galaxies 
that were accreted and disrupted early on, while surviv- 
ing satellites are typically accreted much later (result- 
ing in more metal-enrichment via Type la supernovae) 
and are on average less massive (and experienced less 
efficie nt star formation) than the average halo building 
block (jRobertson et al.li2005t [Font et al.ll2006aD . 

Minor mergers caught in the process of disruption 
bridge the gap between the surviving and completely 
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disru pted dwarf popula tions. The giant southern stream 
(GSS: IIbata et al.ll2001f) in the Andromeda Galaxy (M31) 
offers a rare glimpse at a recent, significant accretion 
event in an external galaxy that is close enough (d ivr.-^i = 
783 kpc; iStanek fc Garnavichl ll998: Holland" iggS") to al- 
low spectroscopy of individual red giant branch (RGB) 
stars. This stream is believed to result from the merger 
with M31 of a g alaxy of stellar ma ss ~ 2 x 10^ Mq less 
than 1 Gyr ago (jFardal et al.ll2007l ). and is estimated to 
have contributed at least 10% of the stars in M31' s halo 
(assum ing azimuth al sym metry of M31's halo; Font et ahl 
l200(nlirbata et al.ll2nn7l ^. 

Ground-based photometric observations have shown 
that the GSS pollutes a large portion of the southern 
quadrant of M31's stellar halo. It is observable to a 
projected distance of i?proj ^ 100 kpc from the cen- 
ter of M31, with a metal-rich "core"^ surround ed by 
a lower metallicity "envelope" (jlbata et al.ll2007f l. The 
linc-of-sight distance to the GSS has been measured from 
^proj ~ 10 to 60 kpc using the apparent magnitude of the 
tip of the RGB. The GSS is located ~ 100 kpc behind 
M31 at i?proj ^ 60 kpc, but is at the same dista nce as 
M31 at i?proj ^ 10 kpc (jMcConnachie e t al.|[200l . The 
detailed star formation history of a field located on the 
GSS has been derived from ultra-deep HST observations 
that reach below the main sequence turn-off; the mean 
age of the stellar population is 8.8 Gyr (70% of the stars 
are younger than 10 Gyr), an d the mean metallicity is 
-0.7 dex (jBrown et al.lf20"06"ah . 

Spectroscopic observations of RGB stars in M31's 
GSS have been limited to the metal-rich core. Line- 
of-sight velocity distributions have been published for 
fields ranging from i?proj = 12 to 60 kpc in projected 
radial distanc e from M31 ' s cen ter (Ibata et al.' 1200^ : 
iGuhathakurta et all 120061 : IKalirai et al. ,2006b). At 
-Rproj = 60 kpc, the GSS's mean line-of-sight (heho- 
centric) velocity is Uhei = —320 km s~^. The inner- 
most unambiguous detection of the GSS is in a field at 



21 kpc, with a mean line-of-sight velocity that 
is blucshifted substa ntially relative to the 60 kpc field: 
Uhoi= -513 km s-i ()Kalirai et al.l I2006bl) . The 12 kpc 
field does not have a secure GSS detection; only 2-4 
stars are found that may be relat ed to the GSS based on 
their velocities (jlbata et al.ll2004 ). In each field in which 
it has been observed, the GSS has a narrow velocity dis- 
persion of ~ 10 to 20 km s~^. 

These existing photometric and spectroscopic data sets 
have inspired detailed models of the collision that pro- 
duced M31's GSS. It is believed to have formed fol- 
lowing a close pericentric passage of an M 31 satellite 



on a hi ghly radial orbit (llbata et a l. 2004: Fon t et aTl 
2006bl: iGeehan et al.l 120061: iFardafet a l. 2006, 12007 
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models showed that several other observed substructures 
in M31 [the Northeast (NE) and Western (W) shelves] 
can be explained as the forward continuation of the GSS, 
and predict a third SE shelf feature that has sinc e been 
discovered spectroscopically |Gilbert et al.ll2007D . The 
latest models demonstrate that several of the observed 
properties of the GSS and other debris in M31, includ- 

In this context "core" refers to the high surface brightness 
portion of the GSS extending along its length, not the GSS's pro- 
genitor, which has yet to be identified. 



ing t he asymmetric surface brig htness distribution of the 
GSS (jMcConnachie et al.ll2003f) , the metal-r ich core and 
metal-poor envelope structure (llbata et"aI1 l2007f) . and 
perhaps even the minor- axis arcs (Ib ata et al.ll2007l ). can 
be explained by a rotating progenitor with a reasonable 
initial metallicity gradient (FaxdiiLet al. 2008). Despite 
these recent successes, the current location of the pro- 
genitor's core and whether or not it is still intact remain 
unknown. More observations of the GSS and related de- 
bris are needed to further constrain models of the colli- 
sion and produce more precise predictions for the current 
location and dynamical state of the progenitor's core. 

As part of the Spectroscopic and Photometric Land- 
scape of Andromeda's Stellar Halo (SPLASH) survey, 
we have conducted an extensive spectroscopic survey 
of M31 RGB stars using t he DEIMQS spectrop:raph 
on the Keck II te l escope ( Guhathakurta et al. 2005]; 
IKalirai et all l2006al : iGilbert et al. .2006. 2007), in clud- 
ing observations of the GSS (jGuhathakurta et al.l [20061 : 
IKalirai et al.l [20b6bD . In this contribution, we present 
spectroscopic results for four new fields (and reanalysis 
of two older fields) in the vicinity of the GSS. We identify 
multiple kinematically cold features in their velocity dis- 
tributions that we associate with the GSS or other photo- 
metric overdensities in M31's stellar halo. These include 
the innermost unambiguous kinematical detection of the 
GSS and the first kinematical detection of GSS "enve- 
lope" debris. We also present the kinematical detection 
of substructure in the region of the photometrically iden- 
tified 60 kpc minor-axis arc (Jbata et al. .2007) , which 
some mod els of the col l ision suggest could be related to 
the GSS (jFardal et al.ll2008D . In addition, we further 
characterize the properties of a secondary ki n ematic ally 
cold component discovered by IKalirai et al.l (j2006b( ) in 
their study of a GSS core field at i?proj = 21 kpc. The 
kinematical and chemical properties of stars in the GSS 
presented here will provide improved observational con- 
straints on the orbit and internal structure of the dwarf 
satellite progenitor and M31's potential. 

The observations, reduction procedures, and our tech- 
nique for identifying a clean sample of M31 RGB stars 
are described in §|2] The properties of stars in individ- 
ual spectroscopic fields are discussed in §|3J the kine- 
matical components in each field are identified and their 
properties (mean velocity and velocity dispersion) are 
measured. The properties of the GSS and other debris 
features are discussed and compared in §|4l Finally, the 
observations are discussed in the broader context of stel- 
lar halos in §|S] and the main results of the paper are 
summarized in §[51 

2. DATA 

In this section, we describe the data reduction and 
analysis techniques that produce the final M31 RGB 
sample analyzed in the following sections. We describe 
the selection of fields from our spectroscopic survey 
(ii l2.1[) . give an overview of our photometric and spectro- 
scopic observations and reductions (§§ I2.2I - [2?3|) . and dis- 
cuss the [Fe/H] measurements used in our survey (S I2.4p . 
Finally, we describe the technique used to isolate the 
general M31 field population of RGB stars from MW 
dwarf stars and M31 dwarf spheroidal members f? l2.5p . 
and present the resulting velocity distributions of stars 

(S. 
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Fig. 1. — Locations (in M31-centric coordinates) of the 
Keck/DEIMOS spectroscopic fields disc ussed in this paper (§[2]|, 
overlaid on a star count map of M31 lllbata et al.| [2b05X The 
GSS is the overdensity of stars extending to the south-southeast 
of M31's center. The approximate size and position angle of each 
DEIMOS spectroscopic slitmask is indicated by a rectangle. Black 
rectangles denote the location of masks from our survey of M31 
for which velocity distributions have been previously published (a3 
and H13s; Guhathakurta ct al. 2006; Kalirai ct al. 2006b), while 
green rectangles denote masks whose velocity distributions are be- 
ing presented for the fi rst time. Red rect angles denote masks IFl, 
IF2, IF6 and IF8 from Ibata et al.' f200l). The polygons show the 
regions chosen by Ibata ct al. (2007) for analysis of the photomet- 
ric properties of the GSS (red solid outline), GSS core (red dashed 
outline), GSS envelope (red dashed-dotted outline), and the 60 kpc 
southeast minor-axis arc "stream C" (blue solid outline). 

2.1. Field Selection 

The data presented here represent a smah subset 
of the fuU data set collected as part of our SPLASH 
Keck/DEIMOS spectroscopic survey of M31 RGB stars. 
The subset was selected based on location of the fields 
with respect to M31's GSS and other debris that may be 
related to the GSS (Fig.[T]). In relatively high-density re- 
gions of M31, "fields" consist of one spectroscopic mask 
(field f207) or two overlapping spectroscopic masks (field 
HI 3s), while in lower density regions fields consist of mul- 
tiple masks designed from a single KPNO 4 m/Mosaic 
pointing (e.g., fields al3 and m4). Fields dl and d3 were 
initially targeted to probe the stellar kinematics of the 
dwarf spheroidals And I and And III (Kalirai et al., in 
preparation) . We are using them here to study the stellar 
kinematics of the underlying M31 field halo population, 
including possible GSS stars. 

Fields f207, H13s, and a3 are located on/near the 
highest-surface brightness region of the GSS core (near 
its eastern edge) at -Rproj ~ 17, 21, and 33 kpc, respec- 
tively. Fields dl and al3 are located to the west of the 
highest-surface brightness region of the GSS at i?proj ~ 
45 and 58 kpc, respectively, but are s till well within th e 
broader region spanned by the GSS (|Ibata et al.ll2007D . 
while field d3, at i?proj ~ 69 kpc, is on the western edge of 



this region. Finally, field m4 is located on M31's south- 
east minor axis at i?proj ~ 57 kpc, overlapping one of the 
minor-axis arc features ( "stream C" ) iden tified as stellar 
streams in the photometric survey of Ibat a et aTl ([2007) . 
The latest models show that for certain collision geome- 
tries, similar arc morphologies on the minor-axis can be 
produced, and thus stream C could be related to the GSS 
(iFardal et an[200l) . 

2.2. Photometric and Spectroscopic Observations 

Photometric catalogs for fields f207 and H13s were 
derived from images taken with the MegaCam instru- 
ment^ on the CFH T 3.6-m telescope (for details, see 
iKalirai et ai]l2006bf ). Ima ges were obtained with the 
g' and i' filters. The observed stellar magnitudes were 
transformed to Johnson-Cousins V and / magnitudes us- 
ing observations of L andolt photometric standard stars 
([Kalirai et al.l [2006af ). Photometry for the remainder 
of the fields (a3, al3, m4, dl, and d3) was derived 
from images taken with the Mosaic Camera on the 
KPNO 4-m telescope^. The KPNO images were ob- 
tained with the Washington system M an d T2 filters 
and the intermediate-width DD051 filter (Osth einierl 
|2003() . The observed magnitudes were transformed to 
V and / magnitudes usin g the transformation equations 
of lMaiewski eTall (poM ). 

The DD051 filter measures the surface-gravity sensi- 
tive Mgb and MgH stellar absorption features. There- 
fore photometry in the intermediate- width DD051 band, 
when combined with photometry in the M and T2 bands, 
can be used as a discriminant of stellar surface-gravity. 
This allows photometric pre-selection of spectroscopic 
targets that are likely to be red giant branch (RGB) stars 
in M31 (as opposed to foreground MW dwarf stars). For 
masks designed from Mosaic photometry, the spectro- 
scopic target selection algorithm gave highest priority to 
objects with high photometric RGB probability [as deter- 
mined from the star's position in the (M— DD051) ver- 
sus (M — T2) color-col or diagram (e.g. JPalma et al.ll2003l : 
iMaiewski et al.l [2005f) ] . star-like morphology (based on 
the DAOPHOT parameters chi and sharp), and with 
Iq magnitudes within ^ 1 mag of the tip of M31's RGB. 
The available space remaining on the mask was then 
packed with lower priority "filler" targets. All spec- 
troscopic targets were restricted to the brightness range 
20 < /o < 22.5. In fields imaged with CFHT/MegaCam, 
spectroscopic targets were selected based on their /q 
magnitude and the SExtractor morphological criterion 
stellarity. Details of the spectroscopic target selec- 
tion procedure for field s with and without D D051 pho- 
tometry are given b y iGuhathakurta et "all (j2006) and 
IKalir ai et al.l ()2006bD , respectively. 

Spectra were obtained using the DEIMOS spec- 
trograph on the Keck II 10-m telescope, using the 
1200 line mm~^ grating, which has a dispersion of 

* MegaPrime/MegaCam is a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope which is 
operated by the National Research Council of Canada, the Insti- 
tut National des Science de I'Univers of the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii. 

^ Kitt Peak National Observatory of the National Optical As- 
tronomy Observatory is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 1 

Details of Spectroscopic Observations and Basic Results. 



Mask 


Projected 


Pointing 


center: 


PA 


Date of 


Photometric 


# Sci. 


# 


# of M31 




Radius 


19000 


<5 T9000 


(°E of N) 


Obs. (UT) 


Bandpasses'' 


targets'" 


Stars''''^ 


Starst'-d 




(k-Dc) 


/h .m .s\ 

v • ■ ) 


io.i.n\ 
\ ' ' } 














f207.1 


17.4 


00:43:42.64 


+40:00:31.6 


-27.0 


2005 Nov 4 




202 


123 


114 


H13s.l= 


21.3 


00:44:14.76 


+39:44:18.2 


+21.0 


2004 Sep 20 




134 


100 


87 


H13s_2= 


21.3 


00:44:14.76 


+39:44:18.2 


-21.0 


2004 Sep 20 




138 


121 


98 


a3_l* 


33.9 


00:48:21.16 


+39:02:39.2 


+64.2 


2002 Aug 16 


M, T2, DD051 


85 


20 


17 


a3 2^ 


32 5 




o y . . ou . o 


+178 2 


9009 Orf 1 1 


i V J J -'^ 2 1 J—' "J ^ 


81 


32 


28 


a3_3^ 


31.9 


00:48:23.1 / 


+oy:12:oo.o 


+270.0 


zUUo Uct ZD 


Ti /T n~' TM^Ot^l 

M , J2i -UJJUoi 


83 


43 


30 


ai_i 


44. 1 


00:45:48.61 


+oo:Uo:4d.o 


0.0 


zOOo INOV 


M , J2i -UJJUoi 


155 


91 


21 


dl 2 


45.6 


nn-4R-1 
uu.iu.-LO.yo 


T^OO.UU.^ i .u 


+90.0 


zuuu oep ±u 


i VJ ; -'^ Z 1 J—' J-' U -1- 


149 


96 


28 


al3.1 


58.5 


00:42:58.34 


+36:59:19.3 


0.0 


2003 Sep 30 


M, r2, DD051 


80 


25 


10 


al3.2 


60.6 


00:41:28.27 


+36:50:19.2 


0.0 


2003 Sep 30 


M, r2, DD051 


71 


23 


10 


al3.3 


56.4 


00:42:25.76 


+37:08:25.4 


0.0 


2005 Nov 5 


M, T2, DD051 


78 


29 


13 


al3.4 


56.5 


00:41:32.15 


+37:08:37.4 


0.0 


2005 Nov 5 


M, r2, DD051 


86 


32 


11 


d3.1 


68.3 


00:36:03.83 


+36:27:27.4 


+90.0 


2005 Sep 8 


M, r2, DD051 


120 


89 


5 


d3.2 


69.8 


00:35:39.61 


+36:21:41.8 


0.0 


2005 Sep 8 


M, T2, DD051 


122 


84 


8 


m4_l 


59.3 


01:00:40.24 


+38:42:07.6 


-90.0 


2006 Sep 20 


M, r2, DD051 


72 


30 


13 


m4_2 


58.3 


01:00:46.90 


+38:51:05.9 


-90.0 


2006 Sep 21 


M, T2, DD051 


81 


22 


5 


m4_3 


57.1 


01:00:50.72 


+39:01:00.6 


-90.0 


2006 Nov 20 


M, T2, DD051 


90 


40 


15 


m4_4 


56.3 


00:59:15.92 


+38:42:12.0 


-88.3 


2006 Nov 20 


M, T2, DD051 


70 


31 


6 


m4_5 


54.9 


00:59:18.67 


+38:52:22.5 


-84.0 


2006 Nov 23 


M, T2, DD051 


75 


38 
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^ The catalogs used for spectroscopic target selection were derived from photometry obtained in the listed photometric bandpasscs f ^ |2.2| l . 
^ A number of targets were observed on two different masks. The total number of unique science targets/stars/M31 RGB stars in fields H13s, 
a3, dl, and d3 is therefore less than the reported number. The number of stars with velocities recovered in two separate observations is 9 in 
H13s, 7 in a3, 3 in dl, and 11 in d3. Of these, 7 stars in H13s, 6 in a3, and 1 in dl are M31 RGB field stars. Of the duplicates in the dwarf 
satellite fields, 1 each in dl and d3 are MW dwarf stars; the remainder are dSph members. 
The number of stellar spectra with secure velocity measurements. 

The number of M31 RGB stars is defined as the number of stars that are identified as secure or marg inal M31 RGB stars by the lGilbert et al.| 
||2006[) diagnostic method (§[23J. For the dl and d3 masks, additional selection criteria have been applied to exclude members of the dSphs 

(§[2XTJ. ^ 

^ Previously published in iKalirai et akl ( |2006bD , the data presented here have undergone new reductions with our improved pipeline (§[23J. 
The sample of M31 RGB stars presented here is 70% larger than in IKalirai et akl ( |2006bD . 

^ Previously published in Guhathaku rta et al.| j2006l) , the data presented here have undergone new reductions with our improved pipeline 
(^ 12.31 1 . The data was originally published without the benefit of the lGilbert et al.| ( 120061) diagnostic method for determining M31 membership 
(§[23J. 

0.33 A pix . The observations were obtained over 
five observing seasons (Fall 2002-2006) and are sum- 
marized in Table [1] The spectra were reduced using 
modified versions of the spec2d and spec Id software de- 
veloped at the University of California, Berkeley in sup- 
port of the DEEP2 galaxy redshift survey^°. These rou- 
tines perform standard spectral reductions: flat-fielding, 
night-sky emission line removal, extraction of the two- 
dimensional spectra, and cross-correlation of the one- 
dimensional spectra with template stellar spectra. The 
template library includes stellar spectra obtained with 
DEIMOS and galaxy templates from the Sloan Digi- 
tal Sky Survey. The cross-correlation procedure has 
been tested extensively a gainst radial velocity stan- 
dards (jSimon fc Gehal[2007l ) . and successfully reproduces 
known velocities. All spectra and cross-correlations were 
visually inspected, and only spectra with secure veloc- 
ity measurements were included in further analysis. The 
observed velocity was corrected both for the motion of 
the earth around the sun (heliocentric correction) and 
for imperfect centering of the star within the slit (based 
on the observed position of the atmospheric A-band 



absorption feat ure relative to ni ght s ky emission lines; 
ISimon & rxehall2007l : ISohn et all [20071 ). The typical er- 
ror in the velocity measurements are ~ 6 km s~^ for 



See http : //astron . berkele y . edu/~coop6r/de6p/spec2d/primer , 
|http : //astron!Ter^eTeyTe3u7~c ooper/deep/specld/primer .h tmT] 



the M31 RGB stars presented here. It is calculated by 
adding in quadrature the random velocity measurement 
error (estimated from the cross-correlation routine) and 
a systematic error of 2.2 km s~^ (estimat ed via repeat 
observations of stars: IsTmon fc Gehall2007f) . 

2.3. Previously Published Fields 

The fields H13s and a 3 ha ve been previously an alyzed 
bv lKalirai et~aLl (|2006b[) and lGuhathakurta et al] ([2006), 
respectively. These were among the first fields observed 
and reduced in our M31 spectroscopic survey. Since the 
original reduction of these fields, we have optimized our 
pipeline for reducing and measuring the velocities of stel- 
lar objects. It now includes a greater number of spec- 
tral templates (obtained with the same instrument as 
the science spectra) for measuring velocities via cross- 
correlation as well as a correction for imperfect centering 
of the star in the slit (ij l2.2|) . These improved data reduc- 
tion and cross-correlation routines result in smaller sys- 
tematic and random velocity errors. Velocities obtained 
with the improved cross-correlation procedure are shifted 
by ^ +20 km s^"'^ from the velocities obtained with the 
original cross-correlation pipeline. The typical measured 
velocity error of M31 RGB stars in these fields has been 
reduced from ~ 15 km s~^ to ~ 6 km s~^. 
_The data in fields a3 and II13s presented here have 



^fffij pn re-reduced using the new cross-correlation routine. 
All papers published by our group from 2007 onwards 
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utihze our improved reduction pipehne. We have also 
recovered a significant fraction of spectra in field H13s 
that had failed in our original reduction pipeline. The 
sample of M31 RGB stars in fi eld H13s is now 70% l arger 
than the sample presented bv IKalirai et al.l (|2006b( ). 

In addition to the improvement in the reduction 
pipeline, the data in field a3 have benefited from an im- 
proved method to separat e M31 RGB stars f rom fore- 
ground MW dwarf stars (|Gilbert et al.l 120061 see ^ [231 
below). This method had not yet been fully developed 
when these data were first presented. It is particularly 
important for securely identifying RGB stars in the kine- 
matically broad spheroid of M31, whose radial velocities 
can overlap the tail of the MW dwarf star velocity dis- 
tribution. 

2.4. Photometric and Spectroscopic Metallicity 
Estimates 

We estimate metallicities using both photometric and 
spectroscop ic techniques f or each star in our M31 sur- 
vey. Like IKalirai et al.l ()2006aD . we estimate photo- 
metric metallicities, [Fe/H]phot, by comparing the star's 
position in the [/q, {V — /)o] color-magnitude diagram 
(CMD) to a finely spaced g rid of theoretical is o chron es 
[t = 12 Gyr, [a/Fe]=0.0; IVandenBerglFan (pOM )]. 
The observed stellar magnitudes are transformed into 
Johnson- Cousins V and / magnitudes using the trans- 
formations described in ^ \2.2l a nd exti nction/reddening 
corrections based on the Schlcgc fet al.l ([1998) dust maps 
are applied on a star-by-star basis. A distance modulus 
of (to — M)o = 24.47 (corresponding to a distance to 
M31 of 783 kpc; iStanek fc GarnavichI Il998f l is applied 
to the isochrones. Although several systematic differ- 
ences in our photometric metallicity estimates could ex- 
ist between the data from the two different filter systems 
pCalirai et al. 2 006a. ) , tests of the data show that any 
such differences are small compared to other sources of 
systematic error in our photometric metallicity estimates 
(discussed below). 

The spectroscopic metallicity estimates, [Fe/Hjspoc, are 
based on the strength of the Ca ii-triplet absorption lines 
at A - 8500A usi ng the empirical calibration relations 
of iRutledge et all (|1997al lbl) and assuming a horizontal 
branch magnitude of Iq = 25.17 for M31 stars. The 
calibration is derived from RGB stars in MW globular 
clusters that span a range of metallicities. 

It is reassuring that in general our photometric metal- 
licity estimates are in good agreement with the spec- 
troscopic metallicity estimates based on the strength of 
the Call absorption line triplet. The individual Can 
EW measurements a nd derived [Fe/H]sr>B r estimates have 
large random errors (jGilbert et al.il20Cl6| ) due to the low 
typical signal-to- noise of our M31 RGB spectra {S/N ^ 
15 per A). Since the [Fe/H]spec estimates are very noisy, 
we can only make a meaningful comparison if we bin our 
sample of stars in terms of [Fc/Hjphot and compute the 
average [Fe/H]spoc in each bin (Fig. 7 of IKalirai et al.l 
[2006a). We thus use photometric rather than spectro- 
scopic metallicity estimates for the analysis presented 
here. 

Although the random error is smaller for [Fe/H]phot 
than [Fe/H]spec, there are substantial systematic errors 
in [Fe/H]phot due to differences between the age and 



alpha-enrichment of a star and the adopted isochrones. 
Moreover, there are substantial variations between the 
isochrones computed by different groups. Photo metric 
metallicity estimates based on the Padova (Girard i et al.l 
I2002D and Yale-Yonsei (Y^; Dcmarq ue et al.i l2004f) 
isochrone set s vary from those based on the Victoria- 
Regina (VR; IVandenBerg et all 2006) isochrones by a 
median amount of ~ 0.1 dex. The variation across the 
isochrone sets is largest for the bluest stars 0.3- 
0.5 dex). 

There are no observational constraints on the alpha- 
enrichment of M31 halo stars, but the assumption of a 
single age for the stellar population is known to be in- 
correct. Deep HST/ACS observations of fields in M31's 
stellar halo and the CSS that reach below the main- 
sequence turnoff have shown that the age of the stel- 
lar population ranges from ~ 8 - 13 Gyr in the halo and 
~ 6-13 Gyr in the CSS; ^ 33% of the population is 
younger than 10 Gyr in M31's halo, while ~ 70% of 
the popula tion is younger than 10 Gyr in the CSS field 
(jBrown et al. 2006a, 2007). These studies find that in 
general the intermediate-age stars are relatively metal- 
rich (-0.75 < [Fe/H] < 0.4) while the older populations 
are relatively metal-poor ([Fe/H] < —1.0). 

Our photometric metallicity estimates are therefore 
likely to be biased low for the most metal-rich stars by 
our assumption of an old (12 Gyr), coeval stellar popula- 
tion, since the HST studies show that these metal-rich 
stars are typically of intermediate-age. Changing the 
isochrone age from 12 to 6 Gyr results in a ~ -1-0.15 dex 
shift in the derived stellar metallicity for stars with 
[Fe/H]phot > —1.0, and a ^ -1-0.3 dex shift for stars with 
—2.3 < [Fe/H]phot < —1.5. Although the systematic bias 
in the derived metallicity estimates introduced by assum- 
ing an old stellar population is large for intermediate-age, 
metal-poor stars, the HST studies show th at such stars 
are rare in both M31's halo and the GSS (jBrown et al.l 
[2006a). 

2.5. Isolation of the Field MSI RGB Sample 

Even in fields with DD051-based photometric prese- 
lection of likely M31 RGB stars (§[221), MW dwarf stars 
along the line of sight to M31 still constitute a signifi- 
cant fraction of the recovered stellar spectra. Since the 
radial velocity distributions of MW dwarf stars and M31 
RGB stars overlap, we use a set of photometric and spec- 
troscopic diagnostics to determine the probability that 
an individual star is a red giant at the distance of M31 
or a dwarf star in the MW. The diagnostics are (1) ra- 
dial velocity, (2) photometry in the Washington M and 
T2 bands and the DD051 band, (3) Nai EW versus 
(y-/)o color, (4) position in an [Iq, {V - I)o] CMD and 
(5) spectroscopic versus photometric [Fe/H] estimates. 
This technique for separating the M31 RGB and MW 
d warf star populations is presented and tested in detail 
bvl Gilbert etall ((2006h . We briefly summarize it here. 

Each of the five diagnostics provides separation be- 
tween M31 RGB stars and MW dwarf stars based on 
different physical parameters. The combination of the 
five diagnostics provides a powerful test of whether an 
individual star is an M31 red giant or MW dwarf. While 
M31 RGB stars have a broad velocity distribution around 
M31's systemic velocity of Vsys = —300 km s~^, MW 
dwarf stars tend to share the sun's rotation about the 
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Fig. 2. — Overall likelihood dist ributions for stars in each of the 
fields analyzed in this paper fS lZSj l. The shaded histo grams denote 
stars identified as M31 RGB stars. Stars with (L^) > that are not 
shaded are stars that, despite having physical characteristics that 
are consistent with the M31 RGB population, have been removed 
from the M31 RGB sample due to their very blue (V — /)o colors. 
In the dl and dS fields, only stars not as sociate d with the dSph 
galaxies are included in this distribution (i) l2.5.1ll . 

MW's center and have Uhoiio ~ km s~^. As discussed 
in §[2ll the combination of the M, T2, and DD051 pho- 
tometry is sensitive to the difference in surface gravity 
between RGB and dwarf stars. The strength of the Nai 
doublet at A 8190A is both surface-gravity and tem- 
perature sensitive. It provides good separation of RGB 
and dwarf stars with (y — /)o > 2. M31 RGB stars follow 
the RGB isochrones in the CMD, while MW dwarf stars 
are more broadly distributed in CMD space. Finally, due 
to the fact that both the [Fe/H]phot and [Fe/H]spoc esti- 
mates are based on the assumption that the star is an 
M31 RGB star, these estimates agree only if the star is 
in fact an RGB star at the distance of M31. 

A star i is given a likelihood in a given diagnos- 
tic j of being an M31 RGB star, defined as Lij = 
log(Pgiant/-Pdwarf ), whcrc P is the probability the star 
is an RGB or dwarf star based on a comparison of the 
star's position in that diagnostic space to probability dis- 
tribution functions based on empirical tr aining sets of 
definite M31 RGB and MW dwarf stars () Gilbert et all 
120061 §3). The individual likelihoods are combined to 
form an overall likelihood {Li). If (i^) > 0, it is more 
probable that the star is an RGB star than a dwarf star, 
and vice versa for {Li) < 0. Figure [2] shows the overall 
likelihood distributions of stars in each of our fields. For 
the purposes of the current analysis, we use stars with 
(Li) > — i.e., more likely to be an M31 RGB star than 
an MW dwarf star. Furthermore, we include only stars 
whose positions in the [Iq, {V — I)q] CMD are within the 
isochrone grid ([Fe/H]phot = —2.31 to 0.0) to within the 
average photometric error of the sample^ ^. 

It is possible that stars with (Li) > that are signifi- 

11 In the terminology of [Gilbert eFall S200B ). we include class 
-1-1, +2, and +3 stars in the following analysis, and exclude class 
—2 stars from our M31 RGB sample. 



cantly bluer than the isochrone grid in the CMD may be 
asymptotic giant branch stars at the distance of M31 or 
M31 RGB stars with flawed photometry. However, MW 
dwarf stars along the line of sight with large negative ve- 
locities {v\os S —150) tend to have blue (T^ — /)o colors, as 
these stars are generally main-sequence turn-off stars in 
the MW's halo. Since several of the diagnostics have lim- 
ited differentiating power between RGB and dwarf stars 
for very blue stars [{V — /)o < -1-2], the (Li) values of 
these stars are strongly influenced by their M31 RGB- 
fike velocities (Gilbert et al. 2006, §4.1.2). We therefore 
choose to be conservative and exclude these very blue 
stars from our M31 RGB sample. These stars represent 
only 2% of the total stel lar sample. 

In lGilbert eFall (pool , we used a large sample of RGB 
stars in M31's inner spheroid (i?proj S 30 kpc) to inves- 
tigate the effect of removing the velocity diagnostic from 
our selection of M31 RGB and MW dwarf stars. This 
sample was ideal for this analysis since the density of 
RGB stars in M31's inner regions is much higher than the 
density of MW stars along the line-of-sight. We found 
that while the exclusion of the velocity diagnostic does 
increase the completeness of the M31 RGB sample at 
velocities between —100 and km s~^, it also increases 
the level of MW dwarf star contamination at these ve- 
locities. Some of the fields presented here are at large 
distances from M31's center (i?proj ~ 60 kpc), where the 
density of MW dwarf stars along the line-of-sight is sig- 
nificantly higher than the density of M31 RGB stars. We 
thus choose lower contamination over higher complete- 
ness for the present work, and include the velocity diag- 
nostic in our selection of M31 RGB stars. This affects 
only stars with velocities in the range of the MW dwarf 
distribution. Since we do not fit the velocity distribu- 
tion of M31's spheroid in this work and the kinematically 
cold substructures presented below have mean velocities 
< -250 km s-i (§[3]), incompleteness of the M31 RGB 
sample at velocities near zero will only affect the esti- 
mated fraction of the stars in M31's kinematically broad 
spheroid. Furthermore, this effect will be small since the 
portion of velocity-space where we expect incompleteness 
due to the use of the velocity diagnostic is both in the 
tail of the M31 spheroid velocity distribution and is small 
compared to the full velocity space used in the fits to the 
velocity distributions in §[3l 

2.5.1. Andromeda Satellite Fields 

The dSph And I is located within the area spanned by 
the GSS, while the dSph And III is located just outside 
this region (Fig. [H marked by the "dl" and "d3" spec- 
troscopic masks, respectively). In addition to dSph RGB 
stars and MW dwarf stars, M31 "field" stars (RGB stars 
that are physically associated with the halo of M31 rather 
than with the dSph) are expected to be present in our 
spectroscopic sample. The half-light radii of the dSphs 
And I and And III are 2.8^ (600 pc) and 1.6' (360 pc), re- 
spectively (iMcConnachie fc Irwinll2"006[ ). therefore much 
of the area of each 4' x 16' DEIMOS slitmask primarily 
targets the M31 field star population. 

M31 field RGB stars are expected to have different 
velocity, metallicity, and spatial distributions than the 
dSph members. Figures a) and[51[a) display the ve- 
locity distributions of all stars observed in the dl and 
d3 fields, respectively (thin open histograms). The bold 
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Fig. 3 . — M31 field star selection in the dl spectroscopic field 
(S I2.5.1|I . (a) Line-of-sight velocity distribution of all stars (thin 
open histogram), and stars iden tified as RGB stars by the likeli- 
hood method described in §[23] (thick, shaded histogram). In ad- 
dition to the prominent kinematically-cold peak formed by dSph 
members, there are RGB stars present in this field with a wide 
range of velocities, as expected for stars associated with MSl's 
spheroid, (b) Sky positions of the RGB stars in field dl. The center 
(red star), half-light radius (inner ellipse) and tidal radius (outer el- 
lipse) of the And I dSph are marked ( McConnachie & Irwin 2006). 
RGB stars at or beyond the tidal radius of the dSph are designated 
as M31 field stars (open diamonds). Open circles denote red giants 
designated as M31 field stars by virtue of their separation from the 
locus of dSph members in the f-[Fe/H]pi,ot plane (panel c). Filled 
circles denote RGB stars that do not meet our selection criteria for 
M31 field stars. The vast majority of these stars are dSph mem- 
bers. Crosses denote MW dwarf stars, (c) Velocity vs. [Fe/H]pijot 
for all stars. RGB stars associated with the dSph form a tight lo- 
cus in this space. A second population of RGB stars with higher 
[Fe/H]pijot values ( [Fe/H] pi^^t > —1.0) and a wider spread in veloc- 
ities, most of which are beyond the tidal radius of the dSph, are 
more likely associated with M31's spheroid, (d) Location of M31 
field stars and dSph members in the [Iq, {V — I)o] CMD. Theoret- 
ical isochrones shifted to the distance of M31 have been overlaid 
on the data. The isochrones are for an age of 1 2 Gyr, [a/Fe] = 0, 
and s pan a range of [Fo/H] from -2.31 to 0.0 I IVandenBerg et alj 
[200^) . 

shaded histograms denote the stars identified as red gi- 
ants using the method described in § 12.51 Although the 
majority of the RGB stars in the two dSph fields are 
part of the prominent kinematically-cold velocity distri- 
butions of the And I and And III dSphs (at ~ —375 
and —345 km s~^, respectively), there are also RGB stars 
spanning a wide range of velocities in each field, as ex- 
pected for stars that are associated with M31's spheroid. 

We utilize the spatial, kinematic, and chemical proper- 
ties of the RGB stars to identify individual RGB stars as 
likely M31 field stars. The spatial distributions of stars 
observed in the dl and d3 fields are displayed in Fig- 
ures 6) and[lj6). The center of each dSph is marked 
by an open star. The inner and outer ellipses denote the 
half-light radii and tidal radii, respectively, of each dSph 
(jMcConnachic & Irwin 2006). All RGB stars beyond the 
tidal radius of the dSph are designated as M31 field stars 
(open diamonds) . It is possible that some of these stars 
were tidally stripped from the dSph, in which case they 
are recent additions to M31's spheroid. 

Figures [3]^ c-d) andUc-d) demonstrate that stars as- 



FlG. 4. — Same as Figure \3\ for field d3 located on the dSph 
And III. 



sociated with the dSph galaxies (filled circles) form fairly 
tight loci in the f-[Fc/H]phot plane and the [/o, {V — I)o] 
CMD. A small fraction of the RGB stars within the tidal 
radii of the dSphs are designated as likely M31 RGB 
stars by virtue of having velocities and/or [Fe/H]phot es- 
timates that are well removed from the dSph locus (open 
circles). All but one of these stars are beyond the half- 
hght radius of the dSph [Figs. M^b) and Wib)]. In the 
d3 field, these stars are removed from the mean velocity 
of the dSph by more than 8 times the measured width 
of the dSph velocity distribution {(Ty). The RGB star 
in the dl field at u = —438 km s~^ that is within the 
tidal radius of the dSph is 5. Say removed from the mean 
velocity of And I. The full RGB population in the dl 
field appears to be comprised primarily of two distinct 
populations: a relatively metal-poor ([Fe/H]phot< —1-0) 
population spanning a narrow velocity range (the And I 
dSph), and a relatively metal-rich ([Fe/II]phot> —1-0) 
population with a broad range of velocities. The stars 
beyond the tidal radius of And I comprise the major- 
ity of the more metal-rich population [Fig. [31[c-rf)], in- 
dicating that this population is not associated with the 
dSph. We therefore designate all RGB stars in dl with 
[Fe/H]phot> -1-0 as likely M31 field stars. 

The selection criteria for M31 field stars described 
above yields an incomplete sample, as M31 field stars 
that fall both within the tidal radius of the dSph and 
near the region of i;-[Fe/H]phot space occupied by the 
dSph members are not included in our sample. However, 
the fraction of M31 field stars that satisfy both these 
criteria is expected to be small, since (i) dSph members 
span a very narrow range in velocity and (ii) the surface 
brightness of each dSph becomes increasingly dominant 
over MSl's spheroid as one moves towards the center of 
the dSphs, resulting in an increasingly greater likelihood 
that spectroscopic slits are placed on a dSph member 
rather than an M31 field star. Only the M31 field stars 
will be analyzed in the current work. The properties of 
the dSph members in these fields will be discussed in a 
future paper (Kalirai et al., in preparation). 
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Fig. 5. — Line-of-sight, heliocentric velocity distributions for stars 
in each of the fields shown in Figure [l] displayed in order of in- 
creasing distance along the giant southern stream (GSS). Th e ve- 
locity distributions of stars identified as M31 RGB stars (§[23]l 
are shown by the bold histograms. The velocity distributions of 
all observed stars (including stars identified as MW dwarfs) are 
shown as thin histograms. Stars associated with the dSphs And I 
and And III have been removed from the velocity distributions. 
Arrows mark the mean velocity of the GSS: black arrows denote 
velocities measured in this work (§[3]l, while grey arrows denote 
previously published velocities. The mean velocity of the GSS be- 
comes increasingly negative as the stream approaches the center 
of M31. In addition to the prominent kinematical peaks identified 
as the GSS, several of our spectroscopic fields show evidence of 
other kinematical substructures (fields f207, H13s, and m4; §[3ll. 
Each panel is labeled with the projected radial distance of the 
field from MSl's center and the total number of stars with radial 
velocity measurements. The majority of the data come from the 
SPLASH Keck/DEIMOS survey of MSl's stellar halo (fields f207, 
H13s, aS, dl, al3, m4, and d3). T he velocities in fields IF8, IF6, 
IF2, and IFl were obtained from llbata et al.l (2004); only stars 
with V < —100 km s~^ are available from the literature (dashed 
lines). M31's systemic velocity is Vaya = —300 km s~^. 



2.6. Velocity Distributions 

Figure [5] presents the line-of-sight hehocentric velocity 
distributions of all stars (M31 red giants and MW dwarfs) 
observed in our fields, as well as the velocity distributions 
of stars in the fou r GSS spectroscopic fields published by 
llbata et all (|200l (fields IFl, IF2, IF6, and IF8). In the 
innermost fields, M31 RGB stars greatly outnumber MW 
dwarf stars along the line-of-sight, while the majority of 
stars observed in the outermost fields are MW dwarf stars 
due to the decreasing density of M31's halo with i?proj- 

Most of the fields show clear evidence of kincmatically 
cold features in their stellar velocity distributions in addi- 
tion to stars with a wide range of velocities (MSl's kinc- 
matically broad spheroid). In fact, several fields appear 
to have multiple kincmatically cold components (e.g., 
f207, H13s, and m4). The GSS is a prominent structure 
in the stellar velocity distributions of fields near the east- 
ern edge of the high surface brightness GSS core (f207, 
H13s, IF6, a3, IF2, and IFl; Fig. [I]). The detection of 
the GSS in field f207 is the innermost secure kinematical 
detection of the GSS to date. There are also prominent 
kincmatically cold peaks in fields dl and al3 with veloc- 
ities in keeping with the expected velocity of the GSS. 



We will show below that these are also attributable to 
the GSS, and are the first kinematical detections of the 
GSS in fields to the west of its prominent, eastern edge 

(§m]& §iiiD. 

3. IDENTIFICATION OF SUBSTRUCTURE IN 
INDIVIDUAL FIELDS 

This section discusses in detail the stellar velocity dis- 
tributions and evidence for the presence of kincmatically 
cold components in each of our Keck/DEIMOS spectro- 
scopic fields (Fig. The velocity and metallicity distri- 
butions of M31 RGB stars are analyzed in order to iden- 
tify substructure in each field, and maximum-likelihood 
fits to the velocity distributions are performed to char- 
acterize the kinematical properties of the substructure, 
including its mean velocity, velocity dispersion, and the 
fraction of stars belonging to it. The fields near the east- 
ern edge of the GSS are discussed first (f207, II13s, and 
a3; ^ I3.2p . followed by the fields to the west (dl, al3, and 
d3; i) l3.3.ip and the field on the 60 kpc minor-axis arc 
(m4; i) l3.3.2p . A broader analysis of the kinematical and 
chemical properties of the identified substructures in re- 
lation to each other and a discussion of the origins of the 
features will be presented in §|31 

3.1. Characterizing Kinematical Substructure 

We use a combination of kinematics and metallicities to 
identify the number of kinematically-cold substructures 
in each field. We then fit a multiple-component Gaus- 
sian to the kinematical distribution using a maximum- 
likelihood technique. The maximum-likelihood multi- 
Gaussian fits in each field are discussed below and sum- 
marized in Table [2] Quoted errors are 90% confidence 
limits from the maximum- likelihood analysis. Although 
the true shapes of the velocity distributions of the vari- 
ous kinematical components in M31's spheroid are likely 
different from pure Gaussians, they provide a convenient 
means of characterizing the mean velocity and velocity 
dispersions of kinematical features. 

We assume a kincmatically broad distribution of ve- 
locities corresponding to M31's well- mixed spheroid is 
present in each field. To date, only two spectroscopic 
studies h ave measured t he ve locity dispersion of M31's 
spheroid. [Gilbert et~aLl ()2007[) used a sample of ~ 1000 
M31 RGB stars in spectroscopic fields from i?proj = 9- 
30 kpc along M31's minor axis. They measured a velocity 
dispersion of o-jpi^ = 129 km s-i for M31's spheroid, and 
saw no evidence for a change in the velocity dispersion 
with radius between 9-30 kpc. The stars were selected 
using the technique des cribed in ^ | 2.5[ but with out use of 
the velocity diagnostic. iChapman et al.l (|2006f ) measured 
an average velocity dispersion of a^P^ = 126 km s~^ 
for M 31's spheroid, consistent with the [Gilbert et al.l 
(j2007l ) value. Their measurement was based on a sam- 
ple of ~ 800 M31 stars selected using windows in line-of- 
sight velocity space designed to remove foreground Milky 
Way stars, M31 disk stars, and substructure. They also 
fit a model in which the velocity dispersion of M31's 
spheroid decreases with increasing projected distance. 
Using this model, they find that M31's spheroid is ex- 
pected to have a velocity dispersion of ct^p'^ = 137 km s~^ 
at i?pj.oj= 17 kpc, decreasing to ct^p'^ = 100 km s^^ at 
-Rproj= 60 kpc. 
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Fig. 6. — Line-of-sight, heliocentric veloci ties of stars in the GSS 
core field at Hproj = 17 kpc (field f207; S lS^TTI l. (a) The bold 
histogram den otes the distribution of stars classified as M31 red 
giants ( ii l2.5l l. The red curve is the maximum-likelihood multi- 
Gaussian fit to the M31 RGB stars (Tabled. The thin histogram 
denotes the velocity distribution of all stars (M31 RGB and MW 
dwarfs), (b) Velocity versus metallicity {[Fe/H]piiot) of the stars in 
field f207. Stars classified as M31 RGB stars are denoted by circles 
and stars classified as MW dwarf stars are denoted by triangles. 

(c) Velocity distribution of M31 RGB stars with [Fe/If]phot > -1.0. 

(d) Velocity distribution of M31 RGB stars with [Fe/If]phot < -1.0. 
The GSS core is evident as a kinematically-cold peak of relatively 
metal-rich stars at d ~ —525 km s" -"^ in panels (a) -(c), and a 
second kinematically-cold, relatively metal-rich group of stars is 
seen at ti ~ —425 km s~^. In addition, stars associated with 
M31's spheroid are present with a wide range of velocities centered 
on the systemic velocity of M31 {vsy. 



-300 km s-i). 



Since the numbe r of M31 spheroid star s at large radii 
(> 40 kpc) in the (iChapman et al.ll2006D study is small 
and a significant portion of velocity space was not used 
for measuring the change in the velocity dispersion with 
radius, we assume a constant value for the dispersion of 
M31's spheroid in analyzing the velocity distribution of 
our fields {{v)''^^^ = -300 km s-\ ct^p^ = 129 km s^i). 
However, for the fields at larger radii (i?proj> 40 kpc), 
we will comment on the effect on the fit of assuming a 
smaller velocity dispersion for M31's spheroid. 

3.2. GSS Core Fields 

3.2.1. New GSS Field at R 17 kpc 

Field f207 is located on the high-surface brightness 
core of the GSS at a projected radius of 17 kpc from 
the center of M31 (Fig. [T|). The hne-of-sight veloc- 
ity distribution of stars in this field is shown in Fig- 
ure [6^ a), while panels (b-d) show the velocity distri- 
bution of stars in field f207 as a function of [Fe/H]phot 
(§ 12. 4p . The GSS is the kinematically-cold, relatively 
metal-rich group of stars at very negative velocities {v < 
—480 km s~^). This is consistent with the expected ve- 
locity of the GSS at this position based on previous kine- 
matical studies (Fig. [5l Hbata et aP 12004 iKalirai et al.l 
l2006bl:[Guhathakurtaet al.ll2006f) . The data in this field 
provide the innermost secure kinematical detection of the 
GSS to date. 

In addition to the GSS, a second prominent kinemati- 



cally cold, metal-rich population can be seen in the veloc- 
ity range —460 < u < —400 km s^^. A less prominent, 
kinematically cold, primarily metal-rich group of stars 
may be present ai v —300 km s^^. As expected for a 
field at this projected radial distance from M31, the ma- 
jority of the stars not in the ki nematically cold peaks are 
also relatively meta l -rich (e.g.. iDurrell et al.ll2001ll2004l : 
IKalirai et al.l[2006al: iGilbert et al.ll2007l) ! 

The data are consistent (via a reduced-x^ test) with 
being drawn from the maximum-likelihood Gaussian 
with three kinematical components: the GSS core at 



= -524 km s 



,-1 



a second kinematically cold com- 



ponent at V — —426 km s~^, and M31's kinematically 
broad spheroid [Table[2l Fig.[6^a), solid curve]. We there- 
fore adopt the triple-Gaussian fit for this field in the anal- 
ysis presented here. Future observations may increase 
the statistical significance of the v ~ —300 km s^^ com- 
ponent, which currently has an excess of only ~ 6 stars 
(~ 5% of the M31 RGB stars in this field) above the 
number expected from M31's broad spheroid. The stars 
in this possible kinematical peak have a mean velocity of 
V —305 km s^^ and a dispersion of 8 km s~^. 

The measured velocity dispersion of the GSS core and 
the second kinematically cold component are cr^^^ — 
23.2tll km s-i and = 20.8l7^g^ km s'^, respec- 
tively. The median velocity error of M31 RGB stars in 
this field is 6.9 km s~^ (5 12. 2p . Thus the intrinsic veloc- 
ity dispersions of the GSS and second kinematically cold 
component are aj,"*'' = 22.2 and cr™*'' = 19.6 km s^^, 
respectively. We will present evidence in (ii l4.3p that the 
second kinematically cold component may itself be a part 
of the GSS. 

3.2.2. Existing GSS Core Fields at 21 kpc and 33 kpc 

Fields H13s and a3 are located near the eastern-edge 
of the GSS core at 21 and 33 kpc (Fig. [l]). Figures [7] 
and [8] display the line-of-sight velocity distributions of 
stars in these fields. Data in HI 3s and a3 have pre- 
viously been presented b y IKalirai et al.l (|2006b| ) and 
'Guha thakurta et al.l ()2006[ ). respectively. The data pre- 
sented in this work benefit from several improvements 
in our reduction and analysis procedures, including an 
improved stellar template library for determining line- 
of-sight velocities, a correction to the measured veloc- 
ity for mis-centering of the star in the slit, and the re- 
covery of data from several CCD chips on the H13s_2 
mask which had initially failed our automated reduction 
pipeline (ij § l2.2l - [273|l . Our revised estimates presented 
below for the mean velocity of the kinematically cold 
components in fields H13s and a3 are consistent with the 
^ -|-20 km s^^ shift we have observed between our im- 
proved and original reduction pipelines (Fig. [5]). 

The two kinematically cold, metal-rich components 
previously identified in field HI 3s are clearly seen in 
Figure [71 along with a kinematically broad, relatively 
metal- rich spheroid component. The parameters for the 
maximum-likelihood triple Gaussian fit to the velocity 
distribution in this field are listed in Table [21 

Although both the GSS core at w ~ -490 km s"^ 
and the second kinematically cold component {v ^ 
—390 km s^^) in field H13s are primarily metal- rich 
[Fig- IZl^jc)], they are also present in the velocity distri- 
bution of stars with [Fe/H]phot < -1-0 [Fig. Wib,d)]. In 
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TABLE 2 

Results of Maximum Likelihood Fits: Fields with One or More Kinematically Cold Components. 



Field 



Primary cold component" 



^j;^ primary 



primary 



Secondary cold component*^ 



Field Location 



Afprimary/A^total (l,>=-°"dary ^secondary N,,,,,,a,,y / Nt,,,i 



£207 
H13s 
a3 
dl 

al3 
m4 



-524.4"' 



4+7.7 
-7.4 

-489.7t*;^ 
-382.0lii;g 



-300.7 
-355.3 



+ 12.0 
11.6 
fS.l 



23.2^ 



,+7.2 
-5.0 

21.3l|0 

+4.6 
-3.3 
13.3 



16.8; 

30.3j 

32. 2j 
11.4 



+5.2 



,+0.08 
^^-0.07 
48+° °'' 

n 5q+° " 
"■^^-0.12 
7+0.15 
-0.17 



0.72 
0.45 



+0.15 
0.21 
+0.14 
0.13 



-388.3t^;* 



-254.i 



t-5.9 
-4.9 



20 

17.0 



+ 12.9 
-7.6 
6.6 
6.9 



1+8.6 

-4.5 



6.6 



+6.0 
2.9 



0.31"' 



0.15"* 



0.16 



GSS core 



0.271°;°^ GSS core 



GSS core 

west of highest-surface 
brightness region of GSS core 
GSS envelope 
+0.10 gg j^p^ minor-axis 
arc (stream C) 



^ Quoted errors are the 90% confidence limits from the maximum likelihood analysis. All fields are assumed to also have a 
broad underlying distribution of stars with parameters {v)"'^^ = —300 km s~^ and al^'^ — 129 km s~^ (|Gilbert et al.lf2007l ). 
^ In addition to the primary kinematically cold component, fields f207, H13s, dl and m4 contain second kinematically cold 
components. 




Fig. 7. — Line-of-sight, heliocentric velocit ies of stars in the 
GSS core field at ijp^oj = 21 kpc (field H13s; ^ 13:2:21 . Panels and 
symbols are the same as in Figure |6] The primary (GSS core) and 
the secondary kinematically cold components at ~ —490 km 
and V ~ —390 km s~^, respectively, are clearly seen in panels (a) - 
(c), as well as stars related to M31's spheroid with a wide range of 
velocities. The shaded boxes in panel (a) denote stars that have 
physical characteristics that are consistent with the M31 RGB star 
population, but whose very blu e {V — I)o colors cause them to be 
classified as MW dwarf stars (^ I2.5| I. 



fact, the GSS and second kinematically cold component 
comprise an estimated 51% and 19% of the M31 RGB 
population with [Fe/H]phot < —1-0, respectively, similar 
to the respective fractions estimated using the full ve- 
locity distribution in H13s (Tabled]). These percentages 
were estimated using a maximum-likelihood fit to the 
[Fe/H]phot < 1-0 velocity distribution in Fig. [^d), with 
all parameters except the fraction of stars in each com- 
p onent held fix e d at th e values in Tabled 

iKalirai et al.l (|2006H l measured dispersions of av = 
19 km for both kinematically cold components in 
this field, which is in agreement with our new measure- 
ments (Tabled). The median velocity measurement error 
is 5.4 km s^^ in field H13s, yielding intrinsic velocity dis- 



FlG. 8. — Line-of-sight, heliocentric veloci ties o f stars in the 
GSS core field at i?proj = 33 kpc (field a3; ii l3.2.2l l. Panels and 
symbols are the same as in Figure |6] The GSS core is evident as a 
relatively metal-rich, kinematically cold population of RGB stars 
at 1) ~ —440 km s~ in panels (a) -(c). The velocity distribution 
of field a3 also has a significant contribution from M31 spheroid 
stars spanning a range of velocities and metallicities. 

persions of tr™*'' = 20.6 km s~^ and cr™*'' = 16.1 km s~^ 
for the GSS core and second kinematically cold compo- 
nent, respectively. 

The prominent kinematically cold, metal-rich peak in 
the velocity distribution of field aS , reported as a kine- 
matic al detection of the GSS by ' Guhathakurta et "all 
()2006f ). is evident in Figure [5] along with a significant 
kinematically broad spheroid component that is more 
metal-poor, as e xpected for M 31 RGB stars at this dis- 
tance from M31 (jKalirai et al.| [2006al. Table [H lists the 
parameters from the maximum-likelihood double Gaus- 
si an fit to the velocity distrib ution. 

iGuhathakurta et al.l (|2006D measured cr^^^ = 21 ± 
7 km s~^ for the GSS in this field, consistent with 
our new measurement of cr^^^ = 16.8133 ■ 
The median velocity measurement error in this field is 
6.8 km s~^, yielding an intrinsic velocity dispersion of 
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cr™*'- = 15.4 km for the GSS core. 

3.3. Fields Off the GSS Core 
3.3.1. Fields to the West of the GSS Core 

South Quadrant Dwarf Spheroidal Field at _R~ 45 kpc — 
Field dl is located to the west of the highest surface- 
brightness region of the GSS (Fig. [T|), yet is still sohdly 
within the region of M31's spheroid dominated by the 
GSS (as ev i dence d by the deeper photometric study of 
llbata et al.l (j2007f )). Although it was designed to study 
the dwarf spheroidal galaxy And I, we have identified a 
total of 48 M31 field stars in addition to the large number 
of dSph members in this field fii l2.5.ip . 

The line-of-sight velocity distribution of M31 field stars 
in field dl is shown in Figure [9] There is a kinematical 
peak of relatively metal-rich stars at v —375 km s~^, 
and a compact group of 9 stars (~ 19% of the M31 
RGB stars in this field) in the v- [Fe/H]phot plane at d ~ 
—290 km s^^ that are even more metal-rich [Fig. 
The spread in metallicity of the stars a,tv^ —290 km s~^ 
is only 0.1 dex. 

M31's halo has been shown to become increasingly 
metal-poor towards its outskirts. The median [Fe/H]phot 
of all the M31 RGB stars in this field is —0.51 dex, signifi- 
cantly higher than the expected metallicity of M31's halo 
at this projected d istance ([Fc/ H] < —1.0: [ Kalirai et al.l 
I2006at iChapman et al. 2006; Ko ch et al.ll2008f ). The se^ 
lection criteria to separate M31 field stars from dSph 
members only precludes stars with [Fe/H]phot < — 1.0 and 
—425 < V < —375 km s~^ from being included in the 
M31 RGB sample f§ l2.5.ip . and thus is not a significant 
factor in increasing the mean metallicity of stars in the 
field. 

The prominent kinematical peaks at v ~ 
—380 km s^^and v ~ —290 km s~^ combined with the 
relatively high metallicity of the stars is strong evidence 
that we have detected substructure in this field. Fur- 
thermore, the data are strongly inconsistent with being 
drawn from a single broad (cr^ > 100 km s^^) Gaussian 
centered on M31's systemic velocity (with a probability 
P^2 ^ 0.1). The parameters of the maximum-likelihood 
triple-Gaussian fit to the M31 RGB stars in this field 
are listed in Table [2l A maximum- likelihood fit with the 
dispersion of M31's spheroid set to the v a lue e xpected 
for_Rproj~ 45 kpc by the [Chapman et al.l p006h study 
(ii l3.1[i returns parameters that are within the 90% 
confidence limits of the fit in Table [D 

The location of field dl (Figs. [1] & [5]) argues strongly 
for the fact that the GSS is the source of at least some 
of the kinematic substructure seen in this field. The 
V ~ —380 km s~^ kinematically cold component has the 
appropriate mean velocity to be the GSS (Fig. O, and 
we will show in ij |4.1l that the kinematical and chemical 
properties of the v ^ —380 km feature are consistent 
with this interpretation. The nine stars with velocities 
~ —290 km s""'^ and [Fe/Hjphot'^ —0.2 dex are evenly 
spread over the regions of spectroscopic masks in which 
we find M31 RGB field stars. They thus represent a 
kinematically cold tidal stream, not a spatially compact 
group of stars such as an intact cluster. 

South Quadrant MSI Halo Field at R ^ kpc — Field 
al3 is located far to the west of the high-surface bright- 
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Fig. 9. — Line-of-sight, iieliocentric velocities of stars in field dl, 
located in the southern quadrant of M31's halo at Rproj = 45 kpc 
(ii l3.3.1l l. Panels and symbols are the same as in Figure |6] and the 
shaded boxes in panel ( a) denote very blue stars as in Figure [7] 
Only stars that are judged unlik ely to be associated with the dSph 
galaxy And I are shown here (§[3]5TTJ- The maximum- likelihood 
triple-Gaussian fit (Table |2]l is shown in panel (a). 

ness core of the GSS at i?proj = 58 kpc, solidly within 
the the prominent "GSS envelope" photometric struc- 
ture (Fig. [1]). The surface brightness of this field, es- 
timated from star counts of M31 RGB and Milky Way 
dwarf stars, is significantly (> 3(t) brighter than expected 
based on the surface brightness profile of M31's stellar 
halo deri ved from star counts in fi elds from 18 kpc to 
165 kpc (jGuhathakurta et al.ll2005f) . This implies that 
field al3 likely contains a significant amount of substruc- 
ture which boosts its surface brightness above what is 
expected for a relatively smooth halo field at a radial 
distance of 58 kpc. 

The line-of-sight velocity distribution of stars in this 
field is displayed in Figure [TOl Stars are present with 
a wide range of metallicity. The relatively metal-rich 
stars ([Fe/H]phot > —1.0) span a fairly limited range of 
velocities (—370 ^ v < —230 km s^^), while the rela- 
tively metal-poor stars ([Fe/H]phot < —1.0) display both 
a wide range of line-of-sight velocities (—550 v < 
— 150 km s""'^) and a peak of stars with a narrower ve- 
locity dispersion centered at w « —300 km s~^. 

The parameters from the maximum-likelihood double- 
Gaussian fit to the full velocity distribution are listed 
in Table [2l If a smaller velocity dispersion of = 
100 km s~^ is assumed for M31's spheroid, the resulting 
maximum-likelihood parameters lie well within the 90% 
confidence limits shown in Table [2] The kinematically 
cold component is estimated to comprise a large fraction 
(~ 70%) of the stars in this field. If the number of stars 
in M31's smooth spheroid is assumed to be 30% of the 
M31 RGB stars in this field, the estimated spheroid sur- 
face brightness for field al3 is fully consistent with the 
observed surface brightness profile of M31's stellar halo. 

The kinematical data are inconsistent {P^2 ^ 0.1%) 
with being drawn from a relatively broad (cr„ > 
100 km s~^) single Gaussian centered on M31's veloc- 
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Fig. 10. — Line-of-sight, heliocentric velocities of stars in the 
south quadrant M31 halo field at Rproj = 58 kpc (field al3; 5 13. 3.11 1 
with the maximum-likelihood double-Gaussian fit (red curve) over- 
laid. Panels and symbols are the same as in Figure [6] and the 
shaded boxes in panel (a) denote very blue stars as in Figure [7] 
The velocity and [Fe/HJphot distribution of M31 RGB stars in this 
field shows evidence of multiple kinematical components. Rela- 
tively metal-poor M31 spheroid stars with a wide range of ve- 
locities are present, as well as a colder kinematical component 
centered at ?; ~ —300 km s"-*^ [panels {b)-{d)]. The relatively 
metal-rich subset [[Fe/H]phot > —1.0, panel (c)] of the stars with 
V ~ —300 km s~^ have a somewhat broader velocity dispersion 
than the relatively metal-poor subset [[Fe/H]pi^o(^ < —1.0, panel 

id)]- 

ity, but are marginally consistent with the maximum- 
likelihood single-Gaussian fit {ay — 7'i't\i km s~^, 
P^2 ~ 6%). However the distribution of the data in 
the [Fe/H]phot-w plane, the enhanced surface brightness 
observed in this field, and its location on a known photo- 
metric substructure (the GSS envelope) provide convinc- 
ing evidence that there is another RGB population in this 
field in addition to M31's kinematically broad spheroid. 
The field's location in the photometrically prominent 
GSS envelope makes this the most likely source of the 
substructure seen in field al3. 

The maximum-likelihood double- Gaussian fit to all the 
stars in this field is a simplification of the distribution 
of the data. Since one of the Gaussian components is 
held fixed with parameters representing M31's kinemati- 
cally broad, smooth spheroid population, only one Gaus- 
sian component is available to fit the triangular distribu- 
tion seen in Figure [TOT 6). Panels (c) and (d) show quite 
clearly the difference in the velocity distributions of stars 
with [Fe/H]phot> -1-0 and [Fe/H]phot< -1-0 dex. Gaus- 
sian fits to the velocity distribution of the stars in each of 
these panels have been overlaid. We assume a broad kine- 
matical component is present in the [Fe/H]phot< —1.0 
population, but not in the [Fe/H]phot> —1.0 popula- 
tion. This is motivated in part by the velocity distri- 
bution of the stars [there are clearly metal-poor stars 
spread over a very large range in velocity in Fig. Iior rf) 
while a similar distribution is not seen in Fig. llOf c)] 
and in part by the fact that MSl's stellar halo has 
been shown to be largely metal-poor at these large radii 
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I2008D . The velocity dispersion of the kinematical peak of 
stars with [Fe/H]phot> —1.0 is = 45t^Q km s~^, and 
the velocity dispersion of the kinematical peak of stars 
with [Fe/H]phot< -1.0 is = Ib.ltli km s"!. These 
values are both marginally consistent with the dispersion 
found in the fit to the full velocity distribution (Table [2]). 
The mean velocities are fully consistent with each other 
and with the results of the fit to the full velocity distri- 
bution. 

We can only speculate for now on the relationship be- 
tween the relatively metal-poor, narrower kinematical 
peak seen in Figure [TU[ d) and the relatively metal-rich, 
slightly broader kinematical peak in Figure llOf c) . An- 
other pointing in a different region of the GSS envelope 
is necessary to determine if this structure in [Fe/H] -vios 
space persists throughout the GSS envelope region, or is 
specific to the location of field al3. If it persists, the in- 
teresting structure seen in the [Fe/H] phot is likely en- 
demic to the stars associated with the GSS envelope and 
represents one coherent substructure. If only part of this 
structure is seen in another location on the GSS enve- 
lope, (i.e., only the relatively metal-poor, narrower kine- 
matical peak or only the relatively metal-rich, slightly 
broader kinematical peak) it would give credence to the 
interpretation of this structure as two independent kine- 
matical substructures in M31's halo which are superim- 
posed along our line-of-sight. For the present, we use 
the fit to the full sample of stars in this field (Table [2]), 
and treat the substructure as one (complex) kinematical 
component, identified as the first kinematical detection 
of the GSS envelope. 

South Quadrant Dwarf Spheroidal Field at 69 kpc — 
Field d3 is far to the west of the high-surface bright- 
ness GSS core, and on the edge of the region of M31's 
stellar halo which appears contaminated by the GSS in 
starcount maps (Fig. [T]). Although the spectroscopic 
masks were designed to study RGB members of the dwarf 
spheroidal galaxy And III. we have also identified a total 
of 13 M31 field stars (§[IXl|). 

The line-of-sight velocity distribution of the 13 M31 
field stars in field d3 is shown in Figure [TT] The 
M31 RGB stars span a wide range in velocity and 
[Fe/H] phot values, and the velocities are consistent with 
being drawn from either a Gaussian with parameters 
(w)«ph= -300 km s-\ alP^= 129 km s'^ or al^^ = 
100 km s^^ [the dispersi on implied for th e proje cted ra- 
dius of this field by the 'Chap man et al.l (|2006 fl study] . 
The maximum-likelihood single Gaussian in this field has 
parameters (v) = -304.8t^i^;^, ay = 69.4t^^?J km s^^. 
Based on the current limited sample of stars, there is no 
kinematical evidence of substructure in this field. 

3.3.2. South-East Minor Axis Field at R ^ 57 kpc 

Field m4 is located on M31's southeast minor axis at 
^proj~ 57 kpc (Fig. [T]), at the same location as the arc- 
like photometric substructure termed 'stream C'. The 
line-of-sight velocity distribution of stars in this field 
is displayed in Figure 1121 Two prominent kinemati- 
cally cold features are present: a relatively metal-rich 
([Fe/H]phot ^ —1.0) group of stars at w ~ —350 km s~^ 
and a secondary group of stars with somewhat lower 
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Fig. 11. — Line-of-sight, heliocentric velocities of stars in the 
dwarf spheroidal field d3, located in the southern quadrant of M31's 
halo at Rproj = 69 kpc (i) l3.3.1l l. Panels and symbols are the same 
as in Figure |6] and the shaded boxes in panel (a) denote very blue 
stars as in Figure [7] Only stars that are unlikely tobe associ- 
ated with the dSph galaxy And III are shown here ('^ 12.5.11 1. The 
maximum-likelihood single Gaussian fit to the data in this field 
has parameters (v) = —305, ay = 69.4 km s~^ (solid curve). A 
Gaussian with parameters (v) = —300, cr„ = 100 km s~^ is also 
shown (dashed curve). 

metallicities a.t v ^ —260 km s~^. Another spectro- 
scopic study of M31 RGB stars ~ 1° southwest of field 
m4 also found evidence of two kinematically-cold com- 
ponents (jChapman et al.l 120081 the mean velocities and 
dispersions of these co mpo nents are designated by the 
arrows in Fig. [T2)) . In ii l4.2l we will show that both kine- 
matically cold components can be identified with photo- 
metric substructure seen in the region of stream C. 

The parameters of the maximum-likelihood triple- 
Gaussian fit to stars in field m4 are presented in Ta- 
ble m The maximum-likelihood triple Gaussian param- 
eters assuming a smaller velocity dispersion for M31's 
spheroid (cr^ = 100 km s~^) are well within the 90% con- 
fidence limits in Table [51 The median velocity measure- 
ment error in field m4 is 5.4 km s^^. Thus the intrinsic 
velocity dispersions of the primary and secondary kine- 
matically cold components are cr"*'' — 10.0 km s^^ and 
3.8 km s~^, respectively; the secondary kinematically- 
cold component is just resolved. 

3.4. Summary of the Results for Individual Fields 

We have detected and measured the mean velocity and 
velocity dispersion of the GSS in 5 fields: f207, H13s, 
a3, dl, and al3. Of these, field f207 provides the in- 
nermost measurement of the mean velocity and velocity 
dispersion of the GSS at i?proj ~ 17 kpc, and field al3 
provides the first measurement of the kinematical prop- 
erties of the GSS envelope. Four of our 7 fields show not 
one but two kinematically cold components in addition 
to the kinematically hot M31 spheroid population. Field 
1207 has a second kinematically cold component with a 
mean velocity ~ 100 km s~^ more positive than that of 
the GSS (i) l3.2.ip . This component mirrors the secondary 
kinematically cold component in field HI 3s, which also 



Fig. 12. — Line-of-sight, heliocentric velocities of stars in th e 
southeast minor axis field at i?proj = 57 kpc (field m4; ^ I3.3.2|l . 
Panels and symbols are the same as in Figure [6] There is a promi- 
nent, kinematically cold peak of metal-rich stars in this field cen- 
tered at t) ~ —360 km s~^, and a second kinematically cold peak 
of more metal-poor stars a,t v ^ —260 km s~^. The arrows denote 
the mean velocities and dispersions of kinematically cold feat ures 
detected in a nearby spectroscopic field llChapman et al1l2008l) . As 
in Figure [7] shaded boxes in panel (a) denote blue stars for which 
our M31 RGB/MW dwarf classification is less certain ('^ 12.51 1. 

has a mean velocity ^ 100 km s^^ more positive than 
the GSS's. The second kinematically-cold components in 
fields 1207 and H13s likely have the same physical origin 
as each other, and may also be part of the GSS fii l4.3p . 
The velocity distribution of field m4 shows two promi- 
nent, kinematically-cold peaks fii l3.3.2p . They confirm a 
previous finding of two kinematically cold components in 
the photometrically-identified minor-axis arc stream C. 

4. A COMPARISON OF THE KINEMATICAL AND 
CHEMICAL ABUNDANCE PROPERTIES OF 
THE GSS AND OTHER DEBRIS 

In §[3l we identified and characterized kinematically 
cold components in individual fields. In this section, we 
widen our focus and discuss the chemical and kinematical 
properties of these components in relation to each other 
and in the context of physical structures known to ex- 
ist in M31's spheroid. We also examine the evidence for 
various physical origins for the kinematical components 
which are not clearly attributable the the GSS: the kine- 
matically cold components in field m4 fii l4.2p and the 
second kinematically cold components in fields f207 and 
H13s (§[4331). 

4.1. GSS Core and Envelope 

Debris from the GSS extends over a large portion of 
the southern-most quadrant of M31's stellar halo. The 
structure of the GSS includes the high-surface bright- 
ness, high metallicity "core" (( [Fe/H] ) = —0.54) sur- 
rounded by a more diffuse, lower met allicity "envelope" 
(( [Fe/H] ) = -0.71) (jlbata et al.ll2007l ). Numerical simu- 
lations of the collision of the GSS's progenitor with M31 
show that the core/envelope structure is naturally pro- 
duced if the progenit or is a d isk galaxy with an internal 
metallicity gradient ([Fardal et al..,2008.) . 
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Fig. 13. — Mean line-of-sight velocity and observed line-of-sight 
velocity dispersion (not corrected for measurement error) of stars 
associated with the GSS a s a fu nction of radial projected distance 
from the center of M31 f ^ 14.1. 11 1. The values are from maximum- 
likelihood, multiple-component Gaussian fits (§|3]l to the stellar 
velocity distributions of each of the fields, including GSS c ore fields 
fope n circles), the GSS envelope field (filled circle), and I bata et al.l 
l i2004{ ) GSS core fields (open triangles). Error bars denote 90% 
confidence limits from the maximum-likelihood analysis. Weighted 
least-squares linear fits to the GSS core data from fields f207, H13s, 
IF6, a3, IF2, and IFl are also shown. The dotted line in the bottom 
panel shows the mean velocity dispersion of these fields. The mean 
line-of-sight velocity and velocity dispersion of the GSS in fields 
dl and al3 (open circle at 3.25° and solid circle, respectively) are 
consistent with the fits to the core fields along the eastern edge of 
the GSS. 

We have presented kinematical detections of the GSS 
core in fields f207 (^ [HXT]) . H13s and a3 f^ lHX^ . We 
have also detected kinematically cold components in 
fields dl and al3 which we attribute to the GSS core 
and GSS envelope fi) l3.3.ip . respectively. We will show 
below that their respective velocity and chemical prop- 
erties corroborate this interpretation. 

4.1.1. Kinematics 

An inspection of the velocity distributions of fields on 
the GSS clearly shows that its velocity (marked with ar- 
rows in Fig. [5]) becomes increasingly positive with in- 
creasing distance from the center of M31 ((w)gss = 
-527 km s-i in field 1207, while (u}gss = -318 km s'^ 
in field IFl). Figure [T3] displays the mean line-of-sight 
velocity and observed velocity dispersion (Table ^ of 
the GSS as a function of projected radius from M31's 
center. The mean v elocity and velocit y dispersion of 
the GSS core in the llbata et al.l (pOOi ) fields (Fig. [5]) 
were calculated using maximum-likelihood double Gaus- 
sian fits, with G^'^^{v) held fixed with parameters (u) — 
-300 km s-\ CT^P'^ = 129 km s"^ (M^- The error bars 
represent the 90% confidence limits from the maximum- 
likelihood analysis. 

We performed a weighted least-squares linear fit to the 
mean line-of-sight velocity of the GSS as a function of 
i?proj- The result, using the fields on the highest surface 
brightness region of the GSS core, is: 



{v) 



GSS 



[63.0(i?proj/l°) - 594.4] kms" 



(1) 



This fit includes the data from fields f207, H13s, IF6, a3, 
IF2, and IFl (Fig. [T]). These fields are represented by 
open circles and triangles in Figure 1131 The open circle 
at 3.25°, field dl, was not included in the above fit and 
is discussed further below. 

Although the line-of-sight velocity of the orbit of the 
GSS progenitor is expected to be strongly non-linear 
over the range of radn observed, the GSS exists because 
the progenitor and stars belonging to the GSS are on 
slightly different orbits. When this is taken into account 
in modelling the interaction which produced the GSS, 
the orbital models produce a roughly linear trend for the 
mean line-of-sight velocity of the GSS as a function of 
radial distance ove r most of the observed distance range 
(|Fardal et al.ll200"6f ). The linear approximation is an ad- 
equate description of the data over the range of radii for 
which the velocity of the GSS has been measured. 

The dl f ield (j^pmi ^ 45 kp c) is very near the region 
chosen by llbata et all ()2007l ) for analysis of the GSS 
core's stellar population and is consistent with an ex- 
trapolation of this region towards M31's center (Fig. [1]). 
The GSS core velocity at the projected radius of our dl 
field, based on the above linear fit, is expected to be 
{v) ^ —390 km s~^. This is consistent with the velocity 
found for the primary kinematically cold component in 
field dl (Table H ? [5XTD . 

In addition to the fields in the GSS core. Figure [T3l also 
displays the mean line-of-sight velocity and velocity dis- 
persion of the kinematically cold component in the al3 
field (filled circle). The mean velocity of the kinemati- 
cally cold component in field al3 is more positive than 
expected from the linear fit to the core fields. The mean 
velocity of the GSS envelope is expected to be relatively 
close to, but not necessarily the same, as the mean ve- 
locity of the GSS core at the same projected radius (M. 
Fardal, private communication). 

There is not a significant trend with -Rpioj in the line- 
of-sight velocity dispersion of stars in the GSS (Fig. [T31 
bottom panel), although the data hints that the veloc- 
ity dispersion in the outer regions of the GSS (i?proj > 
40 kpc) may be slightly larger on average than the in- 
ner regions (-Rproj < 40 kpc). The mean velocity dis- 
persion of the GSS in the fields along its eastern edge 
is {ay) — 20.8 km s~^. The velocity dispersion mea- 
sured in the more western GSS core field dl is consistent 
with this mean within the 90% confidence limits from 
the maximum-likelihood analysis, while the velocity dis- 
persion of the GSS envelope field al3 is within the 99% 
confidence limits. The weighted least-squares fit to the 
eastern GSS core fields is also shown in Figure [T51 The 
velocity dispersion of the GSS in fields dl and al3 is 
consistent with this linear fit. 

4.1.2. MetalUcity Distributions 

Figure [HI shows the position in an [/q, {V — I)o] CMD 
of all stars within ±2ay of the mean velocity of the pri- 
mary kinematical substructure in fields on and off the 
GSS core. The stars within ±2cr^ of the GSS core in 
fields f207, II13s, and a3 have similar distributions in 
color-magnitude space. The stars within ±2ay of the 
kinematically cold component in field al3, located on 
the GSS envelope, are significantly bluer than the stars 
associated with the GSS core. We have not corrected for 
the varying line-of-sight distance to the GSS in FigurefMl 
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Fig. 14. — Position of stars in the [Iq, {V — I)o] CMD for stars 
within it2(T„ of the mean velocity of the p rimary substructure 
in fields on and off the GSS core (§14.1.211. I sochroncs shown 
are 12 Gyr, [«/Fe]= 0.0 from I VandenBerg et al.. (^006) and are 
the same set shown in Figure |3] The most metal-poor (bluest) 
isochrone represents [Fe/H]= —2.31 and the most metal-rich (red- 
dest) isochrone represents [Fe/H]= 0.0. The CMD distribution of 
GSS stars in the three core GSS fields are similar, while the stars 
in the substructure in al3 and the secondary kinematically cold 
component in m4 are on average more blue. 



The difference in the line-of-sight dis tance to the GSS be- 
tween fields a3 and f207 is ~ 20 kpc (jMcConnachie et al.l 
|2003[) . which corresponds to only a 0.05 dex difference in 
apparent magnitude. 

Differences in the CMD distributions of RGB stars are 
caused by differences in the metallicity, age, or alpha- 
enhancement of the stellar populations. F or sufficiently 
old st ellar populations (older than ^ 4 Gvr: lGallart et aLl 
120051) , the RGB is primarily sensitive to metallicity, with 
age and alpha-enrichment playing secondary roles. Al- 
though the alpha-enrichment of RGB stars in M31's 
spheroid and the GSS is not observationally constrained, 
these stars are known to span a range of ages (ij |2.4p . 
in contrast to the 12 Gyr old coeval population we as- 
sume in measuring [Fe/H]phot- With these caveats, we 
interpret differences in the CMD distribution of stars as 
differences in metallicity in the following analysis. Due 
to the age-metallicity degeneracy of the RGB, the differ- 
ences between the metallicity distributions of features, 
which reflect the observed differences in the CMD distri- 
butions, deserve more weight than the absolute metallic- 
ity estimates of the stars. 

The assumed distance to the stellar population also af- 
fects the photometric metallicity estimates of stars. How- 
ever, the line-of-sight distance to the GSS core in individ- 
ual fields is not known with a great deal of certainty, and 
the line-of-sight distances to the other kinematic features 
identified in §[3] are not yet known. We therefore assume 
M31's distance modulus when deriving the photometric 
metallicity estimates. For the GSS kinematic features, 
we will note the effect of assum ing the distance measured 
bv lMcConnachie et al.l (|2003f) on the [Fe/H]phot value. 

The final source of systematic error in deriving the 
mean metallicity abundance of a kinematic feature is the 



Fig. 15. — Metallicity distribution function (MDF) of stars within 
±2(T„ of the mean velocity of the primary substructure in fields on 
and off the GSS core, in histogram (left) and cumulative {right) 
form, (a) MDF of the GSS core components in fields f207 (solid 
histogram and curve), HlSs (dashed histogra m and curve), and aS 
(dotted histogram and curve), discussed in § 14.1.21 (6) MDF of 
stars associated with the kinematically cold component in the dl 
(dotted histogram and curve) and al3 (solid histogram and curve) 
fields. The combined MDF of the GSS core components i n pane l 
(a) is shown for comparison (dashed histogram and curve, § 14.1.21 1. 
(c) MDF of stars associated with the v ^ —350 km s"-"^ (solid 
histogram and curve) and v ~ —250 km s^^ (dotted his togram 
and curve) kinematically cold components in field m4 (§ 14.21 1. As in 
panel (fe), the MDF of the combined GSS core components is shown 
for comparison (dashed histogram and curve). The GSS core MDFs 
from fields f207, H13s, and a3 are consistent with one other and are 
similar to those of the dl field and the m4 primary kinematically 
cold component, while the MDFs of the GSS envelope in field al3 
and the secondary kinematically cold component in field m4 are 
significantly more metal poor. 

assignment of individual stars to kinematic components. 
Individual stars can be associated with the kinemati- 
cally cold components only statistically (i.e., through 
maximum- likelihood Gaussian fits as in §[31). Therefore, 
the [Fe/H] distributions discussed below include stars as- 
sociated with the kinematically cold components as well 
as stars associated with M31's spheroid that happen to 
have velocities within ±2(Ti, of the mean velocity of the 
kinematically cold components. For the GSS core fields, 
this is a relatively small effect since the GSS core is a 
narrow and dominant feature well offset from M31's sys- 
temic velocity (Fig. [S]) . 

Figure [15] displays the [Fe/H] distributions of the stars 
in Figure [TJ] in differential and cumulative form. The 
GSS core has a very similar [Fe/H] distribution in all 
three GSS core fields [Fig. [TSf a)]. with a strong metal- 
rich peak and an extended metal-poor tail. The median 
metallicity of stars within ±.2ay of the mean velocity 
of the GSS component is ([Fe/H])nicd = —0.55, —0.51, 
and —0.44 in fields 1207, H13s, and a3, respectively; 
([Fe/H])„,ed - -0.49 (([Fe/H])„,ean = -0.66 ±0.04) for 
all three fields combined. 

Accounting for the esti mated distance to the GS S in 
each of the above fields (jMcConnachie et al.l 120031 ) re- 
sults in only a ^ —0.04 dex difference in the derived 
[Fe/H]phot estimates, while changing the assumption of 
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the age of the population from 12 Gyr to 8 Gyr results 
in a ~ +0.1 dex difference in the derived [Fe/H]phot esti- 
mates. Based on the multi-Gaussian fits, between 5% 
and 8% (or ~ 3 to 4 stars) of the red giants within 
±2ct„ of the mean velocity of the GSS core in fields 
f207, H13s, and aS are associated with M31's spheroid. 
Within i?proj ~ 30 k pc, MSl's spheroid is found to be 
relatively metal-rich ( [ Kalirai et aLll2006at [Gilbert et al.l 
120071: IKoch et al.ll2008D . However, if we assume that the 
GSS core is more metal-rich than MBl's spheroid, and 
remove the 3 to 4 most metal-poor stars in each field, 
the resulting median [Fe/H]phot estimate for the GSS in- 
creases by only 0.02 dex. 

The [Fe/H] distribution of stars within ±2cr of the 
mean velocity of the —380 km s~^ kinematically cold 
component in field dl is very similar to the [Fe/H] distri- 
bution of stars associated with the GSS core [Fig.fTSTb)]. 
The median metallicity of the stars associated with the 
GSS in field dl is ~ 0.1 dex lower than the median [Fe/H] 
of the GSS core in fields f207, H13s, and a3. If the es- 
timated distance to the GSS in each of the fields is ac- 
counted for, the median metallicity of the GSS in field 
dl is found to be 0.13 dex more metal-poor than in the 
more eastern GSS core fields 1207, H13s and a3. 

Stars within ±2(T„ of the mean velocity of the kine- 
matically cold component in field al3 are on aver- 
age significantly more metal-poor than the GSS core 
[Fig. [T5f 6)]. The median metallicity of the stars associ- 
ated with the kinematically cold component in field al3 
is ~ 0.7 dex lower than the median [Fe/H] of the GSS 
core (([Fe/H])„ed - -1.18, ([Fe/H])„ean = -1.19±0.12). 
This is qualitatively consistent with results from purely 
photometric data that t he GSS envelope i s more metal- 
poor than the GSS core (jlbata et al.ll2007l ). A two-sided 
Kolmogorov-Smirnov (KS) test returns a probability of 
0.001% that the MDF of the kinematically cold compo- 
nent in field al3 and the MDF of the GSS core were 
drawn from the same distribution. If we assume that the 
GSS envelope follows the same trend of increasing line- 
of-sight distance with -Rpioj as the GSS core and correct 
for this effect, the difference in the metallicity between 
the GSS core in fields f207, H13s, and a3 and the GSS en- 
velope in field al3 becomes slightly greater (~ 0.8 dex). 
The largest systematic effect in this field may come from 
the identification of individual stars as part of the kine- 
matically cold component. The double- Gaussian fit to 
the stars in this field (Table ^ implies that 13.5% of 
the stars (~ 5 stars) within ±2ay of the mean velocity 
of the kinematically cold component are associated with 
the spheroid. If we assume that the 5 lowest metallicity 
stars within ±2ay of the kinematically cold component 
are M31 halo stars, the estimated [Fe/H]phot of the GSS 
envelope debris increases by -1-0.25 dex (([Fe/H]}i„cd = 
-0.92, ([Fe/H]),„ea„ = -1.00 ±0.10). 

4.2. Minor Axis Arc at i?proj ^ 57 kpc 

We identified two kinematically cold components in 
field m4: one at v = —355 km which is rela- 
tively metal- rich, and one at v = —255 km which 
is more metal-poor fii l3.3.21 Fig. [T^ . This field is lo- 
cated squarely on the 60 kpc minor-axis arc "stream C " 
(Fig. [T]), identified as a photometric overdensity in star- 
count maps. The M31 starcount map does show evidence 
of two overdensities in the region of stream C: a more 



prominent overdensity that is relatively metal-rich, and 
a less prominent overd ensity that is relatively metal-poor 
(jChapman et al.ll2008D . Both the relative contributions 
of the kinematially-cold components in field m4 (Tablc[2]) 
and their metallicities (Fig. [121 discussed in detail below) 
imply that the v ~ —355 km s^^ component corresponds 
to the former, and the v ~ —255 km s^^ component cor- 
responds to the latter. 

Our field is 1.1° northeast of the str eam C spectro- 
scopic fields F25 and F26 presented by iChapman et al.l 
()2008() . They also find two kinematically cold compo- 
nents in this region: a relatively metal-rich stream 'Cr' 
[{v) = -349.5 km s'^ and crj,'^*' = 5.1 ± 2.5 km s'^ 
based on 9 stars (quoted errors are Icr)], and a rela- 
tively metal-poor stream 'Cp' ((w) ~ —285.6 km s~^ 
and cr"'*'^ = 4.3l|^'4 ^m s^^ based on 6 stars). Arrows in 
Figure [TW a) mark the mean line-of-sight velocities and 
intrinsic velocity dispersions of streams Cr and Cp in 
the F25/F26 field. Due to the s patial proximity of our 
m4 spectroscopic masks and the IChapman et al.l ()2008f ) 
F25 /F26 masks as well as the proximity of the kinemat- 
ical features in velocity space, it is reasonable to identify 
our V = —355 km component with stream 'Cr' and 
our V = —255 km s"^ component with stream 'Cp.' A 
comparison of the mean velocities of the kinematically 
cold components in the m4 and F25/F26 fields suggests 
that their velocities diverge as one moves from the south- 
west to the northeast along stream C. 

The [Iq, {V — I)o] CMD distributions of stars within 
±2(T.„ of the mean velocities of the v = —355 km s~^ and 

V — —255 km s~^ kinematically cold components in field 
m4 are shown in the bottom panel of Figure [TH It is 
clear that on average the stars in the v = —355 km s~^ 
primary component are redder than the stars in the 

V = —255 kms~^ secondary component. The [Fe/H] dis- 
tributions of these stars (assuming a coeval population of 
12 Gyr) is displayed in FigurefTSFc). The differences seen 
in the (/, V — I) CMD translate to clear differences in 
the metallicity distribution functions of the stars. While 
the [Fe/H] distribution of the v = —355 km primary 
cold component in field m4 is relatively metal-rich and 
similar to the [Fe/H] distribution of the GSS, the MDF 
of the V — —255 km s^^ secondary cold component is 
significantly more metal-poor. 

The median metallicity of the v ~ —355 km s~^ 
kinematically cold component is ([Fe/H])i„ed — —0.56 
(([Fe/H])mcan = -0.79 ± 0.12), which is similar to the 
metallicity measu red for stream C of [Fe/H]= —0.6 
(jlbata et al.l 120071 . measured from the region shown 
in Fig. [J)- Usin g 10 Gyr old Padova isochrones 
(|Girardi et al.l[200l ). [Chapman et all (l200l find an av- 
erage metallicity of ([Fe/H]}avg = — 0.74±0.19 for stream 
Cr. 

Based on the double-Gaussian fit to the velocities in 
this field (Table [2|), 14% of the stars within ±2(7y of 
the V ~ —355 km s~^ kinematically cold component in 
field m4 are likely associated with M31's spheroid (~ 3 
stars). M31's spheroid has been found to be metal-p oor 
at large radii (Kali rai et al.l l2006a'. ' Koch et al.|[2008h . If 
we assume that the lowest metallicity stars are those 
associated with M31's spheroid, the estimated median 
[Fe/H] phot of the kinematically cold component increases 
by only +0.02 dex. 
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The V ~ —255 km s^^ kinematically cold component is 
on average more metal-poor than the v ~ —355 km 
component [Fig. I15f c)]. The median metallicity of stars 
within ±2(71, of the mean velocity is ([Fe/H])mcd = —1.33 
(([Fe/H])^ ,,, = -1.31 ±0.18). U sing 10 Gyr old Padova 
isochrones. lChapman et aTl (|2008 ) find an average metal- 
licity of ([Fe/H])mean = -1-26 ± 0.16 for stream Cp. 

A physical association with the GSS? — Based on an anal- 
ysis of all stars withi n the stream C and GSS core regions 
marked in Figure [H llbata et all (|2007[ ) claimed an asso- 
ciation between the GSS and the 60 kpc minor-axis arc 
was unlikely due to significant differences in the proper- 
ties of their stellar populations. However, if the GSS's 
progenitor was a disk galaxy with a reasonable metallic- 
ity gradient, simulations show that arc- like features along 
M31's southeast minor axis with different stellar popu- 
lations from the GSS can be produced during the pro- 
genitor's disruption p^rdal et al.ll2dd8( ). although they 
are farther from M31's center than the observed minor- 
axis arcs. The simulated arc features have mean line-of- 
sight velocities that are close to the systemic velocity of 
M31 (usys — —300 km s~^). There is also observational 
precedence for variations in metallicity along a tidal fea- 
ture: a strong metallicity gradient ha s been observed i n 
the Milky Way's Sagittarius stream (jChou et al.ll2007f ). 
Therefore, it remains possible that the GSS and some of 
the substructure seen in this region share a common pro- 
genitor. Regardless, the relatively high metallicity of the 

V = —355 km s~^, 'Cr' strea m indicates its progenitor 
was likely fairly ma ssive (e.g.. lLarsonlll974l : lMateolll998l : 
IDekelfc Wo^[2003h . 

Conn ection with the extended cluster EC4 — ^Collin s et al. 
( 2009(1 find that stream Cp from the IChapma n et al. 
( 2008") study has a line of sight velocity and metallicity 
consistent with the extended cluster EC4, which was dis- 
covered as part of the CFHT/Megacam survey of M31 
(jibata et al.l |2007[ ) and had RGB members which were 
targeted on the F25/F26 masks. EC4 has a half-ligh t 
radius of 33.7 pc and My = 6.6 (|Mackev et al.ll2006f ). 
placing it solidly in the gulf between globular clusters 
and dwarf spheroidals in size-luminosity space. One pro- 
posed explanation for such objects is that they are the 
stripped cores of dSphs. It is thus tempting to consider 
the possibility that stream Cp is the tidal stream of a 
stripped dSph, whose core remains intact as the clus- 
ter EC4. However, the measured M/L ratio for EC4 
of 6.7 tl%MQ/ Lq ()Collins et al.|[200a ) makes this un- 
likely. ICollins et al.) ((20091) argue that it is also unlikely 
the stream Cp stars are tidally disrupted members of 
EC4 itself, since they are between 15-74 core radii dis- 
tant from EC4 and HST photometry of EC4 shows no 
evidence of distorted isophotes or significant tidal debris 
()Mackev et al.ll2006l ). Our discovery of the continuation 
of stream Cp in a field 1.1° (15 kpc) to the NE con- 
firms that stream Cp is not formed from tidally stripped 
cluster stars. It is most likely that stream Cp and the 

V = —255 km s~^ stream in field m4 represents the tidal 
stream of an unknown progenitor, from which the ex- 
tended cluster EC4 was also tidally stripped. 

In conclusion, we find that the metallicity distribu- 
tion of stars kinematically associated with the v ~ 
—355 km s^^ cold component and the GSS core are sim- 



ilar (Fig. [T5)l . while the metallicity distribution of stars 
associated with the v ~ —255 km s^^ is significantly 
more metal-poor than the GSS core. The two kinemat- 
ically cold components in field m4 are consistent with 
being detections further along st ream C of the kin emat- 
ical substructures 'Cr' and 'Cp' ()Chapman et al.ir2008) . 
The progenitors of these streams remain unknown. 

4.3. Second Kinematically Cold Component in the 
GSS Core Fields 

4.3.1. Radial Trends 

Each of the two innermost GSS core fields (adjacent 
to each other at -Rproj ~ 17 and 21 kpc) has a second 
kinematically cold component in addition to the primary 
kinematical component identified as the GSS core (Fig.[5l 
iiij lXO - [3X2| . The top panel of Figure [16] displays line- 
of-sight velocity vs. i?proj for the M31 RGB stars in fields 
f207 and H13s. Stars within ±2(7^, of the mean velocity of 
the GSS core and second kinematically cold component 
in each field (Table[2]) are denoted by blue and red points, 
respectively. Stars comprising the main GSS peak and 
the second cold component follow nearly parallel linear 
trends in the line-of-sight velocity vs. i?proj plane. These 
trends are continuous across the boundary between the 
two fields. 

The second cold component, like the GSS, is evenly 
distributed across the two fields. The bottom panel of 
Figure [16] shows the ratio of the number of stars in the 
second kinematically cold component to the number in 
the GSS (-/Vscc/-^GSs) as a function of i?proj- There is no 
apparent trend in the ratio Nscc/Nqss with i?proj- 

A more accurate measure of the velocity dispersions of 
the two kinematically cold components can be obtained 
by accounting for the linear trend of line-of-sight veloc- 
ity with -Rproj (Fig [S])- The line-of-sight velocities of 
individual stars are shifted to a reference frame defined 
by the line-of-sight velocity of the GSS as a function 
of R proj- A ma ximum-likelihood , triple-Gaussian fit to 
the data (§ 13. ip is performed. The best-fit values for 
the observed velocity dispersions of the GSS and sec- 
ond kinematically cold component are 20.1^|'^ km s~^ 
and 14.ltg g km s~^, respectively. These measurements 
correspond to intrinsic velocity dispersions of 19.4 and 
13.0 km s^^, respectively, when the mean radial veloc- 
ity measurement error in each field is taken into account 

(§§[323]-[3T2|). 

There is no evidence in the other GSS core fields of a 
continuation of the second kinematically cold component 
seen in fields f207 and H13s (Fig. [5]). If the linear trend 
seen in Figure [16] is extrapolated, the mean velocities 
expected for the feature in fields IF8 (-Rproj ~ 12 kpc), 
IF6 (-Rproj ^ 26 kpc), and a3 (-Rproj 33 kpc) are —476, 
—344 and —278 km s~^, respectively. The field dl does 
have a second kinematically cold component, which is 
also offset by ^ 100 km from the primary (GSS core) 
kinematically cold component. It is possible that the 
second kinematically cold component in fields f207 and 
H13s is a tidal stream that follows a radial line from fields 
f207 and H13s to field dl (and therefore is not present 
in fields IF6 and a3). However, this interpretation is un- 
likely based on the metallicities of the stars. Stars in the 
secondary kinematically cold component in field dl are 
uniformly metal-rich (Fig. [9]), and are significantly more 
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Fig. 16. — Comparison of stars associated with the GSS core and 
the second kinematically cold compone nt in fi elds £207 and H13s, at 
Rproj ~ 17 and 21 kpc, respectively (">j |4.3.1l l. (a) The line-of-sight 
velocity distribution of M31 RGB stars in fields 1207 and H13s as 
a lunction ol projected radial distance from M31's center. Open 
circles denote stars within ±2a-y of the mean velocity of the GSS 
core in each field, while solid circles denote stars within ±2(t„ of the 
mean velocity of the second kinematically cold component in each 
field. Solid triangles denote the remainder of the M31 RGB stars 
in these fields. The solid line is the linear fit to t he mea n velocity 
of the GSS core as a function of radius (Fig. I13l S I4.1.1|I . offset by 
+100 km s~^ for the second kinematically cold component. The 
dashed and dotted lines are the linear least-squares fits to the GSS 
core and second kinematically cold component data, respectively. 
(h) The ratio of the number of stars in the second kinematically cold 
component and the GSS core (Nscc/NQsg) as a function of radius. 
The ratio was calculated using running bins of width 3'. The error 
bars are based on Poisson statistics. There is no apparent trend 
with radius in A^soc/A^GSS- The dashed line represents a constant 
ratio of A^scc /Nqss = 0.70, which is the ratio of stars in the second 
kinematically cold component to stars in the GSS based on the fits 
in Table [2] 

metal-rich than the second kinematically cold component 
in fields f207 and H13s ([Fe/Hjphof^ —0.2 dex, compared 
to [Fe/H]phot~ -0.5 dex). 

4.3.2. Implications for the Origin of the Second 
Kinematically Cold Component 

The discovery of the second kinematically cold com- 
ponent {(v) = —388 km s~^) in the GSS core field at 
i?Droi ~ 21 k p c, field H13s, was originally reported by 



iKalirai et alj (|2006bD . These authors considered three 
possible origins for the secondary kinematically cold com- 
ponent: (i) a feature related to M31's disk, (ii) a separate 
halo tidal stream superimposed on the GSS, and (iii) a 
wrapped around component of the GSS. In this section, 
we re-evaluate the merits of each of these as the possible 
physical origin of this feature in light of the newly ob- 
tained data in the GSS core field f207 (_R£roj ~ 17 kpc), 
the improved statistics in field H13s (i) l2.3|) . and recent 
theoretical and observational developments. We also 
consider a fourth possible or igin which was not consid- 
ered bv IKalirai et al.l (|2006bf) : a bifurcation in velocity- 
space of the GSS. We will argue that the data in fields 
f207 and H13s favor a physical connection with the GSS 
rather than a separate physical origin for the second kine- 
matically cold component in fields f207 and H13s. 



M31's stellar disk — The disk of M31 has been observed 
to large radial distances from the center of M31. An 
extended smooth stellar distribution has been observed 
out to ~ 40 kpc in the disk plane, and substructures 
with disk-like velocities have been observed con tinuing 
out in radius to 70 kpc fe.g.. Ilbata et al.l 120051 ). Fur- 
thermore, the mean line-of-sight velocity of the second 
kinematically cold component is within ~ 50 km s~^ of 
the expected line-of-sight velocity for M31 disk stars in 
circular orbits at the location of our masks. We there- 
fore consider the possibility that the secondary kine- 
matically cold component in fields f207 a nd HI 3s rep- 
resents a population of smooth disk stars. IKalirai et al.l 
(|2006bf) show that the surface brightness of the secondary 
kinematically cold component in field H13s [i?disk '~ 
75 kpc assuming an unwarpe d disk of inclination 77° 
(jWalterbos fc KennicuttI 119881 )] is more than 25 times 
brighter than would be expected based on an extrapola- 
tion of M31's smooth exponential disk, effectively ruling 
out the smooth disk of M31 as the o rigin of the substruc- 
ture (Fig. 11 of IKalirai et"aLll2006bl ). The surface bright- 
ness of the secondary peak in field f207 (i?disk ^ 59 kpc) 
is equivalent to that of the secondary kinematically cold 
component in field H13s, and is approximately 10 times 
brighter than the extrapolation of M31's smooth disk. 
The surface brightness of M31's disk decreases exponen- 
tially with increasing projected radius while the surface 
brightness of the GSS is relatively constant over the ra- 
dial range covered by the f207 and H13s spectroscopic 
masks. The 1207 and H13s masks combined cover a dis- 
tance of ^ 30 kpc in t he plane of the disk, which i s 
more than 5 scale lengths (jWalterbos fc Kennicuttl[l988l ) . 
Thus, the ratio of the number of stars in the second kine- 
matically cold component to the number of stars in the 
GSS would be expected to decrease substantially with 
increasing projected radius if the second kinematically 
cold component was part of M31's disk. This is contrary 
to what is observed [Fig. [TBf 6)]. These arguments rule 
out the possibility that the origin of the second kinemat- 
ically cold component in these fields is M31's extended 
stellar disk. 

Unrelated substructure — C osmological s i mula- 
tions of halo formation (e.g., iBuUock et al.l 120011 : 
iBuUock fc JohnstonI 120051 ) as well as observational 
studies of the Milky Way (e.g., IMaiewski et al.l 119961 : 
INewberg et al.l 120021: IBelokurov et all 120061 ) suggest 
that multiple tidal streams may exist along a single 
line of sight. The secondary kinematical peaks in fields 
f207 and H13s could represent a tidal stream in M31's 
halo unrelated to the progenitor of the GSS. Also, as 
mentioned above, substructures rotating with disk-like 
velocities hav e been observed as far as ^ 70 kpc from 
M31's center (jibata et al.ll2005D . The arguments against 
a smooth disk origin for the second kinematically cold 
component in fields f207 and H13s do not preclude the 
possibility that it is substructure related to the disk. 
However, for either halo or disk-related substructures 
we would not expect a direct correlation between the 
line-of-sight velocity of the GSS and the second stream. 
Since the two kinematical substructures are superim- 
posed on each other in projected location, only a limited 
set of orbits would produce the observed correspondence 
in mean hne-of-sight velocity with i?proj (Fig. [T6)) . 



MSl's Giant Southern Stream 



19 



Wrapped around GSS component— iKalirai et all ()2006bD 
discussed the possibihty that the second kinematically 
cold component in field H13s represented part of a later 
orbital phase of the GSS's progenitor that has wrapped 
around M31 and back to the location of this field. 
A few of the many dynamical models of the merging 
event (jlbata et al.ll200l [lb^rteraLll2006bt iFardal et all 
|2006[ ) yield a progenitor orbit that passes through the 
region of fie l ds H I 3s and f207 twice (e.g., Fig. 8 of 
Fardal et al.l 1200^. This origin was disfavored by 
Kalirai et al.l (|2006bl ) because the stars are predicted 
to be moving in very different phases during the two 
passages. This results in a predicted line-of-sight ve- 
locity offset between the stellar streams produced by 
each phase that is inconsistent with the observed veloc- 
ity offset between the GSS core and the second kine- 
matically cold component in fields H13s and 1207 (Ta- 
ble [2]) . The most recent dynamical models of the colli- 
sion suggest that the Northeast and Western shelves ob- 
served in M31 are the forward continuation of the GSS 
(|Fardal et al.l 120071 ) . This model is strongly supported 
by recent spectrosco pic and HS T/ACS observati ons of 
the shelves in M31 (iBrown et a l. 2006b; Gilber t et all 
I2007t iRichardson et al.l l2008l )~Although in this model 
the progenitor's orbit does not pass through the region 
ol fields H13s and 1207 a second time, the shelves them- 
selves could be responsible for the second kinematically 
cold component. The extent of the W shelf is obscured 
by the southwest region of M31's disk. It could produce 
a velocity peak with a 100 km s~^ separation from 
the GSS velocity in fields 1207 and H13s, although the 
second velocity peak characteristic ol shell features (one 
peak corresponds to outfiowing and one to inlalling de- 
bris) would have to be depopulated to match the data 
(M. Fardal, private communicaiton) . The W shell would 
have to span almost 180° in position angle to be the 
source ol the kinematically cold component in fields 1207 
and H13s. Future spectroscopic observations will test 
whether the velocity ol W shelf stars are consistent with 
the second kinematically cold component in the GSS core 
fields f207 and H13s. 

Bifurcation of the GSS — The lact that the velocities ol 
stars in the GSS and second kinematically cold compo- 
nent follow the same linear trend with projected distance 
from M31's center, and that the ratio of stars in the 
two kinematical features is approximately constant with 
-Rproj , suggest the possibility of a direct physical associa- 
tion between the two features. The formation of velocity 
caustics during the tidal disruption of an in-falling satel- 
lite can lead to bifurcations in the observed linc-o f-sight 
veloc i ties of a t idal feature (e.g., [Mcrrificld & K uiikenI 
119981: IFardal et al.il2007f ). Thus, the second kinemati- 
cally cold component in fields f207 and H13s could be 
part of the GSS itself (i.e., stars, like those in the GSS, 
that are approaching their first pericentric passage since 
the progenitor's initial encounter with M31). However, it 
may be dilficult to produce a bifurcation with a separa- 
tion of line-of-sight velocities as large as ~ 100 km s~^, 
as seen in fields f207 and H13s. Further modeling of 
the collision will be needed to determine whether veloc- 
ity caustics of this magnitude can be produced in this 
region of the GSS. 

If the second kinematically cold component is related 
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Fig. 17. — Cumulative [Fe/H] distribution of stars within ±2cr„ 
of the mean velocity of the second kinematically cold components 
in fields H13s (dashed curve) and f207 (dotted curve). The cumu- 
lative [Fe/H] distribution of the GSS core (based on the primary 
kinematically cold component in fields f207 and HI 3s and the kine- 
matically cold co mpo nent in the aS field) is shown for comparison 
(solid curve; Fig. I15I I. The second kinematically cold component 
in fields f207 and H13s is as metal-rich as the GSS core. 

to the GSS, either as a bifurcation of the GSS itself or 
the forward continuation of the GSS, the metallicity dis- 
tributions of stars in the GSS and the second kinemati- 
cally cold components should be similar. Figure [T7] dis- 
plays the cumulative [Fe/H] distributions ol stars in the 
GSS core (solid curve) and the second kinematically cold 
component in field 1207 (dotted curve) and H13s (dashed 
curve). Although the second kinematically cold compo- 
nents have a small number of stars more metal-rich than 
any associated with the GSS core, the [Fe/H] distribu- 
tions are not inconsistent. A two-sided KS test returns 
probabihties of 15% and 22% that the [Fe/H] distribu- 
tions of the second kinematically cold components in f207 
and H13s, respectively, were drawn from the same parent 
distribution the GSS core was drawn from. 

In conclusion, the data rule out M31's extended, 
smooth disk as the origin of the second kinematically cold 
components in fields f207 and H13s. While substructure 
in M31 unrelated to the GSS cannot be ruled out by the 
data, the possibility of a connection with the GSS, cither 
via a bifurcation in velocity space of the GSS itself or a 
wrapped around component of the GSS, is an intriguing 
possibility. The common trend of mean line-of-sight ve- 
locity with projected radius of the GSS and the second 
kinematically cold component combined with the similar 
metallicity distributions of their stars would be a remark- 
able coincidence if the substructures are unrelated, but 
are naturally explained if the second component is also 
part of the GSS. 

4.4. Summary of Debris Properties 

We have detected the metal-rich GSS core in fields 
f207, H13s, a3, and dl. Although field dl is located on 
the dSph And I, we are still able to isolate M31 field stars 
from dSph members, and to identify kinematically cold 
features in the M31 field star population. The kinemati- 
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cal and chemical properties of the —380 km s^"'^ kinemat- 
ically cold component in field dl are entirely consistent 
with its identification as GSS core debris. The kinemat- 
ical and chemical properties of the relatively metal-poor 
—300 km s^^ cold component in field af 3 are consistent 
with what is expected for the GSS envelope based on 
empirical evidence from photometric studies and theo- 
retical modeling of the collision that produced the GSS. 
The second kinematically cold component in fields f207 
and Hf 3s has a similar [Fe/H]phot distribution to the GSS 
core, follows the same parallel trend in the linc-of-sight 
velocity vs. i?proj plane, and may also be associated with 
the GSS. 

The kinematically cold fea tures observed in field m4 
confirm the identification by 'Chap man et al.l ()2008[ ) of 
two kinematically cold components (one metal-rich and 
one metal-poor) in the photometric stream C minor-axis 
feature. The origin of the more metal-rich kinematical 
feature is not known. The more metal-poor feature may 
be physically associated with the extended stellar cluster 
EC4, although it is unlikely to be a tidal stream origi- 
nating from EC4 itself. 

5. DISCUSSION 

Mixed stellar populations appear to be an inevitable 
feature of the M31 halo. Many fields in the SPLASH 
spectroscopic survey show conclusive (Table [2]) or sug- 
gestive evidence of multiple components within the M31 
field halo population itself. In addition, fields that target 
M31's dSph galaxies contain dSph member RGB stars as 
well as M31 field halo RGB stars f^ [^XT|) . All fields con- 
tain a non-negligible fraction of foreground Milky Way 
dwarf star and background galaxy contaminants. It is 
useful to keep the notion of mixed populations in mind 
while interpreting star-count maps. In this section we 
discuss several implications of the mixed stellar popu- 
lations discovered in M31 for our understanding of the 
structure and formation of M31's stellar halo. 

Although many of the spectroscopic fields presented 
here specifically targeted tidal debris (f207, H13s, IF8, 
IF6, IF2, and IFl), others were chosen for DEIMOS 
spectroscopy without prior knowledge of the existence of 
tidal debris at their locations (fields a3, m4, al3, dl, and 
d3). Nevertheless, the majority of stars are found to be 
located in distinct clumps in i>ios-[Fe/H] space. More- 
over, many fields contain multiple distinct groupings of 
stars. If the outer regions of M31's stellar halo were 
well-mixed, stars from numerous progenitors would be 
present along any given line of sight, leading to stars that 
are br oadly distributed in uios-[Fe/H] space. iTremaini 
(|1993D calculates that the filling factor of tidal streams, 
/a, in a Milky- Way type halo built from the disruption 
of accreted satellites should be small (/a < 1) at radii of 
~ 50 kpc (a given volume of the halo is dominated by 
stars from one to a few progenitors), while well- mixed 
systems have /s ^ 1 (stars from numerous progenitors 
are present in any given volume). The fields presented 
above indicate that a large fraction of the stars in M31's 
halo at these radii belong to groups of stars that can be 
identified as individua l tidal streams. As discussed by 
iMaiewski et~all (|1996f) . this may prove problematic for 
drawing conclusions about the global properties of M31's 
stellar halo based on limited lines-of-sight. Further stud- 
ies in other quadrants of M31's stellar halo will be able 



to determine whether the dominance of distinct groups 
of stars in these fields is a ubiquitous property of M3rs 
stellar halo, or is mainly a result of the large amounts of 
debris deposited by the progenitor of the GSS. 

The number of distinct kinematical components 
present in many of the M31 stellar halo fields necessi- 
tates careful interpretation of the properties of specific 
regions and tidal features measured from photometric 
studies of M31's halo population. Even if the percent- 
age of the population associated with a tidal feature 
is known via separate kinematical studies, photomet- 
ric studies cannot identify which stars are more likely 
associated with a specific feature rather than the un- 
derlying halo of M31. Metallicity distributions of fea- 
tures derived purely from photometric observations will 
therefore be a sum of the metallicity distribution of 
the feature itself, M31's halo, and any other structural 
(sub)component or feature present in the analyzed re- 
gion. The minor-axis arc stream C provides a prime ex- 
ample of the complexity that can be present in a small 
region of M31's halo (§[3X2] & §113). Furthermore, the 
extensive presence of mixed stellar populations in M31's 
stellar halo, seen in both kinematical (this work : also 
see llbata e t al. 2004, 2005; Guhathaku rta et all 120061 : 
iKalirai etal... 2006b.: .Gilbert et al..,2007.) and photomet- 
ric surveys (e.g., [Ibata et aL n200lt iFereuson et al.l[200l 
llbata et al.l l2007f ). make identifving control M31 halo 
fields for photometric studies problematic. Kinematical 
data provide the ability to identify which stars are likely 
to be associated with an individual (sub)component or 
tidal feature; this provides a distinct advantage over pho- 
tometric studies in the analysis of the properties of indi- 
vidual stellar (sub)components in M31. 

A variety of spatial structures are expected in stellar 
halos that form via the merging of smaller stellar sys- 
tems, including still-bound satellites, spatially coherent 
tidal streams, larger, more diffuse clouds of debris, and 
a spatially smooth, well-mixed stellar population. The 
data presented here display evidence of a correlation be- 
tween kinematical and spatial structure. The M31 dSphs 
And I and And III are very kinematically cold, with low 
velocity dispersions ((T„ < 10 km s~^; Kalirai et al. in 
preparation). As surviving dSphs, these objects are also 
spatially compact. Many of the kinematically cold com- 
ponents observed in our spectroscopic fields (i.e., m4 and 
the GSS core in f207, H13s, and a3) can be identified with 
spatially sharp debris streams. The kinematically cold 
component in field al3 is warmer (has a broader velocity 
distribution) than the kinematically cold components in 
the GSS core fields. It is identified with the GSS enve- 
lope debris, which appears to be a broad cloud of debris 
rather than a sharp tidal stream. Our fields also contain 
a kinematically hot population with a very broad spread 
in velocities; this population represents M31's underly- 
ing, well-mixed stellar halo. 

Metal-rich tidal streams are expected to dominate over 
metal-poor streams in surface brightness limited stel- 
lar halo surveys, due primarily to the fact that more 
massive dwarf galaxies are both more metal- rich (e.g., 
ILarsonlll97l IMateollT998l : IDekel fc W00II2OOI and tend 
to produce hig her sur face brightness stellar streams 
poimston et al . 2008; G ilbert et aLll2009f) . We currently 
have < 50 M31 RGB stars per spectroscopic field in 
M31's outer halo (i?proj ft 35 kpc). A high surface bright- 
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ness tidal stream is easily identifiable in our survey in the 
form of a kinematically cold component seen against the 
dynamically hot underlying stellar halo. Tidal streams 
with a low level of contrast against M31's spheroid are 
difficult to detect due to our limited sample size. We 
therefore expect the kinematically cold components iden- 
tified in our M31 fields to be more prominent among the 
metal-rich than metal-poor stars. In general this is in- 
deed the case, as shown in panels {b-d) of Figs. [6l-fT2l 

6. CONCLUSIONS 

In this study, we presented data for five new spectro- 
scopic fields observed during the course of our SPLASH 
Keck/DEIMOS survey. We also presented improved re- 
ductions and sample selectio n for two of our previously 
pubhshed fields (a3 a nd H13s: [Guhathakurta et al.ll2006l : 
iKalirai et al.l l2006b[ ) and compiled all previously pub- 
lished line-of-sight velocity data in fields on the CSS. 
This combined spectroscopic dataset was used to ana- 
lyze the kinematical and chemical properties of the CSS. 

The main findings of this paper are: 

• We have detected and measured the mean velocity 
and velocity dispersion of the GSS in two new GSS 
core fields: (i) field f207, at i?proj'^ 17 kpc (the 
innermost kinematic measurement of the GSS to 
date), and (ii) field dl at i?proj'^ 45 kpc (the first 
detection of GSS core debris in a field well removed 
from the sharp eastern edge of the GSS). 

• We have made the first kinematical detection and 
measurement of the mean velocity and velocity 
dispersion of the GSS envelope in field al3, at 
-Rpioj^ 58 kpc. 

• We have presented updated measurements of the 
mean velocity and velocity dispersion of two previ- 
ously published GSS core fields at -Rproj~ 21 and 
33 kpc using improved reductions of the spectro- 
scopic data. 

• The full dataset of new and previously published 
line-of-sight velocities in GSS fields was used to 
measure the dependence of the mean line-of-sight 
velocity and velocity dispersion of the GSS as a 
function of projected distance from M31's center. 

• We utilized the line-of-sight velocities in the 
SPLASH survey fields to identify stars belonging 
to the GSS. The metallicity distribution of stars in 
the GSS is found to be similar in all GSS core fields, 
while the mean metallicity of stars in the GSS en- 
velope is significantly more metal-poor (~ 0.7 dex) 
than the GSS core. 



• We have confirmed the presence of two kinemat- 
ically cold components, one relatively metal-rich 
and one relatively metal-poor, in the stream C re- 
gion at -Rproj~ 60 kpc on the minor axis. A com- 
parison of the mean velocities of the kinematic fea- 
tures in field m4 with a previously published field 
1° to the southwest indicates that the observed 
line-of-sight velocities of the two features diverge 
as one moves from the southwest to the northeast 
along stream C. 

• In our innermost GSS core field (-Rpioj'^ 17 kpc) 
we discovered the continuation of a second kine- 
matically cold component previously seen in a GSS 
core field at -Rproj~ 21 kpc. The kinematical and 
chemical properties and number of stars in this sec- 
ond kinematically cold component are tightly cor- 
related with those of the GSS over Ai?pi.oj = 7 kpc, 
indicating that a direct physical association be- 
tween it and the GSS is likely. The second kine- 
matically cold component in these two fields may 
represent a bifurcation in the line-of-sight velocities 
of stars in the GSS. 

Previous spectroscopic data on the GSS h as been pre- 
sentcd piecemeal in several differe nt p apers flbata et ajj 
l2004; GuhathakurtieLall 120061 : iKalirai et al. 2006b). 
The compilation of previously published data and the 
new GSS fields presented here provide a clear portrait of 
the kinematical and chemical properties of the GSS along 
a significant portion of both its length and breadth. This 
data provides ingredients for improved detailed modeling 
of the collision of the progenitor of the GSS with M31, 
and will provide additional observational constraints on 
M31's potential as well as the properties and orbit of the 
progenitor of the GSS. 
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